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High-resolution molecular spectroscopy is a sensitive probe for violations of fundamental symme- 
tries. Symmetry violation searches often require, or are enhanced by, the application of an electric 
' field to the system under investigation. This typically precludes the study of molecular ions due to 

^S) ' their inherent acceleration under these conditions. Circumventing this problem would be of great 

benefit to the high-resolution molecular spectroscopy community since ions allow for simple trap- 
^ ' ping and long interrogation times, two desirable qualities for precision measurements. Our proposed 

solution is to apply an electric field that rotates at radio frequencies. We discuss considerations for 
QQ , experimental design as well eis challenges in performing precision spectroscopic measurements in 

rapidly time-varying electric fields. Ongoing molecular spectroscopy work that could benefit from 
our approach is summarized. In particular, we detail how spectroscopy on a trapped diatomic molec- 
ular ion with a ground or metastable ''Ai level could prove to be a sensitive probe for a permanent 
electron electric dipole moment (eEDM). 
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O ■ I. INTRODUCTION 

• 1— I " 

A. High- Resolution Molecular Spectroscopy as a Probe of Fundamental Physics 



The quest to verify the most basic laws of nature, and then to search for deviations from them, is an ongoing challenge 
at the frontier of precision metrology. To this end, high resolution spectroscopy experiments have made significant 
contributions over the years. For example, the coupling strengths and transition energies between atomic and molecular 
levels are predominantly determined by the electromagnetic interaction. However, the Standard Model does include 
\ fundamental processes, e.g. the weak interaction which have spectroscopic signatures that are both theoretically 
^\ • calculable and experimentally detectable. Parity-violating transition amplitudes, forbidden by the electromagnetic 
' interaction but allowed in the presence of the weak interaction, have been calculated and measured in atomic cesium 0, 
, 0] and ytterbium [J with sufficient precision to test electroweak theory at the ~ 1% level. In addition, high-resolution 
' molecular spectroscopy experiments are underway to probe parity violation in chiral polyatomic molecules [5|-j8i] and 
I to probe nuclear spin-dependent parity violation in diatomic molecules d, Looking outside of the Standard 

Model, precision molecular spectroscopy experiments have been designed to search for time- variation of fundamental 
constants, such as the electron-to-proton mass ratio (Tll - fTsj and the fine structure constant [H!, [3, as well as to 
search for simultaneous parity and time-reversal symmetry violation in the form of permanent electric dipole moments 
(EDMs) 

In most cases, atoms and molecules that are either neutral or ionic can be studied in an effort to observe the same 
underlying physics; however, typically there are technical advantages to selecting one system over the other. Systems 
of neutral, as opposed to ionic, particles are attractive for precision spectroscopic studies due to the relative ease 
of constructing high-flux neutral particle beams, the relatively weak interactions between neutral particles, and the 
lack of coupling between the translational motion of neutral particles and external electromagnetic fields. Conversely, 
charged particles are favored due to the relative ease of constructing ion traps and the long interrogation times 
that come with studying trapped particles. Indeed, some of the most stringent tests of the Standard Model have 
been performed using trapped ions 31-'34'|, and spectroscopy on trapped molecular ions is of fundamental interest 
for studying interstellar chemistry ,35r,37.] . Looking to combine the techniques of ion trapping and high-resolution 
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TABLE I: Theoretical predictions of the electron electric dipole moment, de- Current listings are taken from Ref. [gJ], which 
extracted the numbers from Refs. [stI. Is^. 



CP Violating Model 


de [e cm] 


Standard Model 


\d, \ < 10-^» 


Supersymmetric models 


|de| < W-'^'' 


Left-right symmetric models 


10-2« < |de| < 10-2'5 


Higgs models 


10"^* < |de| < 3 X 10-2'^ 


Lepton flavor changing models 


10-2'' < |de| < 10-2'' 



molecular spectroscopy, several research groups are working to develop experimental platforms for studying ensembles 
of trapped molecular ions t38-43j. 

The additional degrees-of-freedom afforded to molecular systems, in comparison with simple atomic systems, pro- 
vide additional interaction mechanisms and correspondingly more routes for experimental investigation. For example, 
molecular levels are inherently more sensitive to applied electric fields due to the presence of nearby states of opposite 
parity, e.g. rotational levels and/or A-doublet levels. On the surface, this means that the Stark shifts observed in 
molecular spectra will be significantly larger than the corresponding shifts to atomic transitions. More fundamen- 
tally, this means that in relative weak electric fields the quantum eigenstates of an atomic system are still dominated 
by a single parity eigenstate, while the quantum eigenstates of molecular systems asymptotically approach an equal 
admixture of even and odd parity eigenstates. There are several classes of atomic and molecular symmetry violation 
experiments where larger Stark mixing amplitudes give rise to larger signals. For example, the parity violation signals 
already attained in atomic systems [2tJ| are expected to be exceeded by the next-generation of experiments using 
polarized diatomic molecules [loj . Similarly, in experiments designed to search for permanent electric dipole mo- 
ments, the expected signal size scales with the ability to thoroughly mix parity eigenstates and increases dramatically 
when going from atoms to diatomic molecules |4^-|46[. 

Herein lies the conundrum for symmetry violation searches using trapped molecular ions: the electric field required 
to fully polarize the molecules will interfere with the electromagnetic fields necessary for trapping the ions with the 
likely result of accelerating the ions out of the trap. Our solution to this problem is to apply an electric field that 
rotates at radio frequencies. Under these conditions, the ions will still accelerate, however they will undergo circular 
motion similar to charged particles in a Penning trap fsil-fssj or storage ring j47l - [5^ . The nuances of performing 
high-resolution electron spin resonance spectroscopy in this environment will be the main focus of this work, with the 
ultimate goal of demonstrating that such an experiment on the valence electrons in a ground or metastable ^Ai level 
could prove to be a sensitive probe for a permanent electron electric dipole moment (eEDM). 



B. Motivation for Electric Dipole Moment Searches 



The powerful techniques of spin resonance spectroscopy, as applied to electrons, muons, nuclei, and atoms, have 
made possible exquisitely precise measurements of electric and magnetic dipole moments. These measurements in 
turn represent some of the most stringent tests of existing theory, as well as some of the most sensitive probes for new 
particle physics. As an example, the recent improved measurement of the electron's magnetic moment [33| agrees with 
the predictions [s^] of quantum electrodynamics out to four-loop corrections. Compared to the electron work, muonic 
g-2 measurements [U [s^l are less accurate but are nonetheless more sensitive (due to the muon's greater mass) to 
physics beyond the Standard Model. Digging a new-physics signal out of the muon g-2 measurement is made difficult 
by uncertainty in the hadronic contributions to the Standard Model prediction j55l |. One of the primary motivations 
for experimental searches for electric dipole moments (EDM) is the absence of such Standard Model backgrounds to 
complicate the interpretation of these studies. In the case of the electron, for example, the Standard Model predicts 
an electric dipole moment less than 10"^^ e cm [56|. The natural scale of the electron electric dipole moment (eEDM) 
predicted by supersymmetric models is 10~^^ to 10"^^ e cm (57l - [59| (Table H]). The current experimental limit is 
14 1 < 1.6 X 10"^^ e cm feO]. With predictions of new physics separated by nine orders of magnitude from those 
of "old" physics, and with the current experimental situation such that a factor-of-ten improvement in sensitivity 
would carve deeply into the predictions of supersymmetry, an improved measurement of the eEDM is a tempting 
experimental goal. In this paper we will describe an ongoing experiment that we believe will be able to improve on 
the existing experimental upper limit for an eEDM by a factor of thirty in a day of integration time. 
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C. A Brief Overview of the JILA Experiment 



Our JILA eEDM experiment will be based on electron spin resonance (ESR) spectroscopy in a sample of trapped 
diatomic molecular ions. We will use an A-doubled molecular state that can be polarized in the lab frame with a 
lab frame electric field of only a few volts/cm. The very large internal electric field of the molecule, coupled with 
relativistic effects near the nucleus of a heavy atom, will lead to a large effective electric field, £es, on the electron 
spin. Confining the molecules in a trap leads to the possibility of very long coherence times and therefore high 
sensitivity. Trapping of neutral molecules has been experimentally realized recently, but it remains an extremely 
difficult undertaking. Conversely, trapping of molecular ions is straight forward to implement with long-established 
technology. 

On the face of it, measuring the electric dipole moment of a charged object is problematic. Even for a relatively 
polarizable object like a molecule, one must apply sufficient electric field to mix energy eigenstates of opposite parity. 
This field will cause the ion to accelerate in the lab-frame and limit trapping time. We will circumvent this problem 
via the application of a rotating electric bias field, which will drive the ion in a circular orbit. The rotation rate will 
be slow enough that the molecule's polarization can adiabatically follow the electric field, but rapid enough that the 
orbit diameter is small compared to the trap size. The ESR spectroscopy will be performed in the rotating frame. We 
note that this approach is conceptually related to efforts measuring electric dipole moments of charged particles in 
storage rings |47h52| |. but in our case the radius of the circular trajectory will be measured in millimeters, not meters. 
Precision spectroscopy in time-varying fields can be afflicted with novel sources of decoherence and systematic error, 
which will be discussed in Sees. IIV| IVl and IVII 
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D. A Comparative Survey of Ongoing Experimental Worlc 

The primary purpose of this section will be to review experimental searches for eEDM. We will make no attempt to 
survey the rapidly increasing diversity of low-energy [gl] and astrophysical searches for physics beyond the Standard 
Model. A subset of that broad area of endeavor is the search for permanent electric dipole moments (EDMs), and 
a subset within that focuses on electrons (eEDMs). For comparative surveys of the discovery potential of various 
EDM studies see [63l - [67[ . we summarize here by saying that from the point of view of new physics, experiments on 
leptons provide physics constraints complementary to those on diatomic atoms and to those directl y on bare nucleons 
and nuclei. As for the lepton experiments, there is work on the tau lepton [gll, on muons [Sllls^ and of course 
on electrons as discussed in some detail below. The current best neutron EDM measurement was done at ILL 
there are ongoing neutron EDM searches [69l - [7l| . Beam-line measurements on bare nucleons are envisioned 
The current best atomic dipole measurement is an experiment is in the diamagnetic species, Hg, by the Washington 
group [l^. Many other groups are looking for EDMs in diamagnetic (that is, net electron spin S = 0), ground-state 
electronic levels in Hg [73], Xe [Ti - lTTj . Rn [t^, Yb [zi] and Ra |67i, .80-82]. Experiments on diamagnetic atoms 
(with net electron spin 5* = 0) are sensitive to new physics predominantly via the nucleonic contribution to the Schiff 
moment of the corresponding atomic nucleus. Higher-order contributions from eEDM contribute to the atomic EDM 
of S* = atoms [s^, but these are probably too small to provide a competitive eEDM limit. 

For 20 years the most stringent limits on the eEDM have been the atomic-beam experiments of Commins' group 
at Berkeley [gO, HH, I13| ■ That work set a standard against which one can compare ongoing and proposed experiments 
to improve the limit. Here is a brief survey of ongoing experiments of which we are aware. 

For evaluating the sensitivity of an eEDM experiment the key figure-of-merit is ScgT^/N, where £cS is the effective 
electric field on the unpaired electron, r is the coherence time of the resonance, and N is the number of spin-flips 
that can be counted in some reasonable experimental integration time, for instance one week. The statistics-limited 
sensitivity to the eEDM is just the inverse of our figure-of-merit. We will discuss the three terms in order. 

The conceptually simplest version of an eEDM experiment would simply be to measure the spin-fiip frequency of 
a free electron in an electric field Suh, <^d = o^e^iab, where de is the electric dipole moment of the electron [85|. Alas, 
a free electron in a large electric field would not stay still long enough for one to make a careful measurement of its 
spin-fiip frequency; in practice all eEDM experiments involve heavy atoms with unpaired electron spins. An applied 
laboratory electric field distorts the atomic wavefunction, and the eEDM contribution to the atomic spin-fiip frequency 
cod is enhanced by relativistic effects occurring near the high-Z nucleus j86l.l87j. so that ojd = deSeS, where the effective 
electric field £off can be many times larger than the laboratory electric field fiab- The enhancement factor is roughly 
proportional to although details of the atomic structure come into play such that the enhancement factors for 
thallium {Z = 81) and cesium {Z = 55) are —585 ^] and -1-114 |89|], respectively. Practical DC electric fields in 
a laboratory vacuum are limited by electric breakdown to about 10^ V/cm. The Commins experiment used a very 
high-Z atom, thallium, and achieved an £cS of about 7 x 10^ V/cm [60|. There have been proposed a number of 
experiments in cesium ^90.-.92l ] that expect to achieve £eS of about 10 V/cm. A completed experiment at Amherst [ll] 
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achieved £cS = 4.6 x 10^ V/cm in Cs by using £u\, = 4 kV/cm. 

It was pointed out by Sandars [44l - |4a | that much larger S^s can be achieved in polar diatomic molecules. In these 
experiments, the atomic wavefunctions of the high-Z atom are distorted by the effects of a molecular bond, typically 
to a much lighter partner atom, rather than by a laboratory electric field. One still applies a laboratory electric field, 
but it need be only large enough to align the polar molecule in the lab frame. The Imperial College group |22i] is 
working with YbF, for which the asymptotic value of Se S is 2 6 GV/cm [2!. l93 - l99l | . The Yale group [iTHlsil l uses PbO, 
with an asymptotic value of £cS — 25 GV/cm j20|, [2l[ J"l00| . The Oklahoma group has proposed to work with 
PbF, which has a limiting value of £es ~ 29 GV/cni0^. The ACME collaboration ^ will use ThO, with S^g 
~ 100 GV/cm [28[. The Michigan group is working with WC, with £cS — 54 GV/cm We will discus s can didate 
molecules for our experiment in Sec. IIIBI we anticipate having an S^s of around 25 to 90 GV/cm [28l. [sol . [lOl| . 

After feff I the next most important quantity for comparison is the coherence time r, which determines the linewidth 
in the spectroscopic measurement of ujd- In Commins' beams experiment, t was limited by transit ti me t o 2.4 ms. 
Future beams experiments may do better with a longer beam hue [13], or with a decelerated beam [l02| . Groups 
working in laser-cooled cesium anticipate coherence times of around 1 s, using either a fountain or an optical trap 
[91I [9^ . The PbO experiment has r limited to 80 fis by spontaneous decay of the metastable electronic level in which 



they perform their ESR. Coherence in ThO experiment will be limited by the excited-state lifetime to 2 ms [27[. A 
now discontinued experiment at Amherst [93| achieved t = 15 ms in a vapor cell with coated walls and a buffer gas. 
The JILA experiment will work with trapped ions. The mechanisms that will limit the coherence time in our trapped 
ions are discussed in Sees. IIVI andlVl We anticipate a value in the vicinity of 300 ms. 

The quantity £es converts a hypothetical value of de into a frequency uJd, and r sets the experimental linewidth of 
LOd- The final component of the overall figure-of- merit is VN, which, assuming good initial polarization, good final- 
state sensitivity, and low background counts, determines the fractional precision by which we can split the resonance 
line. Since we have defined N as the number of spin flips counted, detection efficiency is already folded into the 
quantity. Vapor-cell experiments such as those at Amherst or Yale can achieve very high values of effective N , atomic 
beams machines are usually somewhat lower, and molecular beams usually lower yet (due to greater multiplicity of 
thermally occupied states.) Atomic fountains and atomic traps have still lower count rates, but the worst performers 
in this category are ion traps. The JILA experiment may trap as few as 100 ions at a time, and observe only 4 
transitions in a second. 

The discussion above is summarized in Table |lll To improve on the experiment of Commins, it is necessary to 
do significantly better in at least one of the three main components of the figure-of- merit. The various ongoing or 
proposed eEDM experiments can be sorted into categories according to the component or components in which they 
represent a potential improvement over the Commins' benchmark. The prospects of large improvements in both r 
and £off put JILA's experiment in its own category. This combination means that our resonance linewidth, expressed 
in units of a potential eEDM shift, will be 10^ times narrower than was Commins'. Splitting our resonance line by 
even a factor of 100 could lead to an improved limit on the eEDM. This is an advantage we absolutely must have, 
because by choosing to work with trapped, charged molecules, we have guaranteed that our count rate, N, will be far 
smaller than those of any of the competing experiments. 

W e note that there are in addition ongoing experiments attempting to measure the eEDM in solid-state systems |l03l - 
Il06l | . These experiments may also realize very high sensitivity, but because they are not strictly speaking spectroscopic 
measurements, it is not easy to compare them to the other proposals by means of the same figure-of-merit. 

Finally, atoms with diamagnetic ground states may have S ^ metastable states amenable to an eEDM search 

[l07l |. Closely spaced opposite parity s tates in Ra can give rise to an S^s on the electron spin larger [108] than in Tl 
or Cs, but very short coherence times jl08j | may make complicate efforts [63| to measure the eEDM in Ra. 



E. Outline 



A brief overview on the molecular level structure where the eEDM will be measured and on how the measurement 
will be performed is given below in Sec.|TTl Some aspects of the experimental design, including production of molecular 
ions and ion trapping will be covered in Sec. IIIIl Difficulties in performing precision spectroscopy in time-varying 
and inhomogeneous electric and magnetic fields will be discussed in Sec. IIVI This will include discussions of trap 
imperfections, stray magnetic fields, and effects of rotating bias fields. Experimental chops used to minimize systematic 
errors will also be explained. In Sec. |Vl the effects on spin coherence time and systematic errors of ion-ion collisions 
will be investigated. An estimate for experimental sensitivity to the eEDM will be given in Sec. I VII The Appendix 
gives a listing of variables used throughout the paper and a sample set of experimental parameters. 
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TABLE II: Figure-of-merit comparison between several recently completed and ongoing eEDM experiments. For ongoing 
experiments these numbers are subject to change and are often order-ol-magnitude estimates. For the JILA entry M is Hf, Th, 
or Pt and x is H or F. 



Group 


Refs. 


Species 


flab [V/cm] 


£cfi [V/cm] 


r [s] 


N [s-i] 


Berkeley 


[60] 


Tl 


1.23 X 10^ 


7 X 10^ 


2.4 X 10-3 


10^ 


Amherst 


[93] 


Cs 


4 X 10^ 


4.6 X 10^ 


1.5 X 10"^ 




LBNL 


[90] 


Cs 


10^ 


10^ 


1 


10^ 


Texas 


[92] 


Cs 


10^ 


10^ 


1 




Penn State 


[9^ 


Cs 


10^ 


10^ 


1 




Yale 


[17-21] 


PbO 


10 


2.5 X 10^° 


8 X 10"^ 




Imperial 


[22, 23] 


YbF 


8.3 X 10^ 


1.3 X 10"' 


10-3 




Oklahoma 


[24-26] 


PbF 


7 X 10^ 


2.9 X 10^'' 






ACME 


[27, 28j 


ThO 


10^ 


10" 


2 X 10-3 


10^ 


Michigan 


[29] 


WC 




5.4 X 10^° 


10-3 




JILA 


This work 


Mx+ 


5 


3 - 9 X 10^° 


0.2 - 1 


~ 10 



II. MOLECULAR STRUCTURE AND THE BASIC SPECTROSCOPIC IDEA 

A. Molecular Notation 

As we prepare this paper, we have not made a final decision as to which molecule we will use. For reasons discussed 
below, the main candidates are diatomic molecular ions Mx"*", where M = Hf, Pt, or Th and x = H or F. In the case 
of molecules such as HfF+, ab initio methods [13, Il0l| enable us to determine that the '^A state is well described 
by a set of Hund's case (a) quantum numbers: J,S, E, A, $7, Mj, e//. Here J is the sum of electronic plus rotational 
angular momentum, S the total electronic spin angular momentum, E the projection of S onto the molecular axis, 
A the projection of L, the electronic orbital angular momentum, onto the molecular axis, and the projection of J 
onto the molecular axis. In a case (a) '^A molecule jfij can take the values one, two or three. Mj is the projection of 
J along the quantization axis and the labels e// specify the parity of the molecular state. 

In addition to these quantum numbers, the experiment will be concerned with the nuclear spin quantum number 
/, the total angular momentum quantum number F, given by the vector sum of J and /, and mp the projection of 
F along the quantization axis. Throughout this paper we shall assume a total nuclear spin of / = 1/2, the nuclear 
spin of fluorine or hydrogen. This leads to the values F = 3/2 and = 1/2 for the states of experimental interest. 

B. Choosing a Molecule 

In selecting a molecular ion for this experiment we have several criteria. First, we want a simple spectrum. Ideally, 
we would like the supersonic expansion to be able to cool the molecules into a single internal quantum state so that 
every trapped molecule could contribute to the contrast of the spectroscopic transition. Failing that, we want to 
minimize the partition function by using a molecule with a large rotational constant, most likely a diatomic molecule 
with one of its atoms being relatively light. Small or vanishing nuclear spin is to be preferred, as are atoms with only 
one abundant isotope. Second, we need to be able to make the molecule. This requirement favors more deeply bound 
molecules and is the main reason we anticipate working with fluorides rather than hydrides. Third, the molecule 
should be polarizable with a small applied electric field, i.e. it should have a relatively small A-doublet splitting, ujoi- 
Fourth, and most important, the molecule should have unpaired electron spin that experiences a large value of £cff ■ 

These latter two requirements would appear to be mutually exclusive: a small A-doublet splitting requires a 
large electronic orbital angular momentum, which prohibits good overlap with the nucleus required for a large £cS- 
Fortunately, working with two valance electrons in a triplet state allows us to satisfy our needs. One valance electron 
can carry a large orbital angular momentum making the molecule easily polarizable, while the other can carry zero 
orbital angular momentum gi ving it good overlap with the nucleus and generating a large feff- This concept was 
detailed by some of us in Ref. (SO) and for the '^Ai state of i nter est here, the two valance electrons occupy molecular 
a and S orbitals. Our calculations, as well as those of Ref. [lOl| . indicate that in the ^Ai state of ThF"*" and HfF+ 
we should expect uj^i < 27r x 40 kHz with £off ~ 90 GV/cm for ThF+ and £cff ~ 30 GV/cm for HfF+ [28l[loH. 
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C. 1^1 = 1 vs. \n\ = 3 

We mention one final valuable feature we look for in a candidate molecule: a small magnetic g-factor, so as to 
reduce the vulnerability to decoherence and systematic errors arising from magnetic fields. To the extent that spin- 
orbit mixing does not mix other |f2| = 1 states into a nominally "^Ai molecular level, it will have a very small magnetic 
moment, a feature shared by PbF in the ^Ili/2 state [l^. This is because E = —A/2, and because the spin g-factor 
is ^2 times the orbital g-factor. Under these conditions, the contributions of the electronic spin and orbital angular 
momentum to the net molecular magnetic dipole moment nominally cancel. In HfF"*", the magnetic moment of a 
stretched magnetic sublevel level of the '^Ai, J = 1 rotational ground state is about 0.05 /is. This is a factor of 20 
less than the magnetic moment of ground state atomic cesium. In the ^A3 level, on the other hand, the magnetic 
moment in the stretched zeeman level is 4.0 ^b- The \ — 3 state may nonetheless be of scientific interest. The '^Ai 
and •^As levels have ScS equal in magnitude but opposite in sign. If one could accurately measure the science signal, ujd, 
in the '^Aa level despite its larger sensitivity to magnetic field background (and despite its shorter spontaneous-decay 
lifetime), the comparison with the ^Ai result would allow one to reject many systematic errors. 

D. |n| = 1, J = 1 A-doublet 

Since we have not made a final decision as to which molecule we will use, and also because we have yet to measure 
the hyperfine constants of our candidate molecules, the discussion of level schemes in this section will be qualitative 
in nature, usually emphasizing general properties shared by all the molecules we are investigating. To simplify the 
discussion, we will specialize to discussing spectroscopy within the J — 1 rotational manifold of a molecular ■^Ai level. 

For Hunds' case (a) molecular levels with |A| > 1, each rotational level is a A-doublet, that is, it consists of two 
closely spaced levels of opposite parity. We can think of the even (odd) parity level as the symmetric (antisymmetric) 
superposition of the electronic angular momentum lying predominantly parallel and antiparallel to the molecular axis 
[Fig.[lja)]. The parity doublet is split by the A-doubling energy cjcf. A polar diatomic molecule will have a permanent 
electric dipole moment, dmf , aligned along the internuclear axis h, but in states of good parity, there will be vanishing 
expectation value (n) in the lab frame. An applied laboratory electric field, Srot, will act on dmf to mix the states 
of good parity. In the limit of dmf^rot 3> ujcf, energy eigenstates will have nonvanishing (n) in the lab frame. More 
to the point, fl, a signed quantity given by the projection of the electron angular momentum on the molecular axis, 
(L + S) -h, can also have a nonzero expectation value [Fig. [T](b)]. Heuristically, it is the large electric fields developed 
internal to the molecule, along n, that gives rise to the large value of £eS that the electron spin can experience in 
polar molecules. In the absence of the A-doublet mechanism for polarizing the molecule, a much larger field would 
be necessary, dmfi^iab ^ 2i?e, to mix rotational states with splitting typically twice the rotational constant B,,. For 
HfF+, we estimate ojcf will be 2tt x 10 kHz, whereas will be about 2tt x 10 GHz. For a dipole moment = 4.3 D, 
mixing the A-doublet levels will take a field well under 1 V/cm, whereas "brute force" mixing of rotational levels 
would require around 10 kV/cm. For an experiment on trapped ions, the smaller electric fields are essential. 

In the context of their eEDM experiment on the a'^Ei level in PbO, DeMille and his colleagues have explored in 
some detail ^17nl9j the convenient features of an |r2| = 1, J — 1 state, especially with respect to the suppression of 
systematic error. Our proposal liberally borrows from those ideas. In a molecule with at least one high-Z atom, ■^Ai 
states will be very similar to the a'^Ei state of PbO, but with typically smaller values of LOd and much smaller values 
of magnetic g-factor. Singly charged molecules with spin triplet states will necessarily have an odd-Z atom, and thus 
the unavoidable complication of hyperfine structure, not present in PbO. 

In Fig. [T] we present the '^Ai, J = 1 state with hyperfine splitting due to the fiuorine 1=1/2 nucleus. A key feature 
is the existence of two near-identical pairs of mi^-levels with opposite parity. As seen in Fig. [TJb), an external electric 
field, Srot, mixes these opposite parity states to yield pairs of mj^-levels with opposite sign of £es [3] relative to 
the external field. Fig. [T](c) shows the effect of a rotating magnetic bias field, parallel with the electric field, applied 
to break a degeneracy as described in Sec. IIVDI below. Note that any two levels connected by arrows in Fig. [TJc) 
transform into each other under time reversal. Time reversal takes mp — >■ —mp, ^ — ^ — and B — —B, where B 
is the magnetic field. If we measure the resonant frequency for the transition indicated by the solid (or dashed) line 
once before and once after inverting the direction of the magnetic field, time reversal invariance tells us the difference 
between the two measurements should be zero. In the presence of an eEDM, which violates time-reversal invariance, 
this energy difference W'^{B) — W'^{—B) will give 2de£cS- As well, under the same magnetic field the transitions 
indicated by the solid and dashed lines should be degenerate, if the magnetic g-factors are identical for the states 
involved |l09l |. With non-zero eEDM the energy difference VF" — W \ also gives 2de£ee- 

Potential additional shifts, due predominantly to Berry's phase |llOj . are discussed in Sec. IIVI but for now we 
note only that in the absence of new physics (such as a nonzero eEDM) the energy levels of a molecule in time- 
varying electromagnetic fields obey time-reversal symmetry. Reversing the direction of the electric field rotation while 
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(a) (b) (c) 




F = 1/2 — 

FIG. 1: Energy levels of HfF"*" in the ^Ai, J — 1 state including hyperfine structure associated with the fluorine 1—1/2 nucleus. 
A and E are defined as the projection along the molecular axis of the electronic orbital angular momentum, and spin, respectively. 

= A + S. (a) In zero electric field, the eigenstates of the system are states of good parity, |e) — = +f) — |r2 = —l))/y/2 
and I/) — {\Q — +1) + |n = — f))/\/2, separated by a small A-doublet splitting, (b) An electric field, Srot, mixes the parity 
eigenstates yielding states with well defined Q. (c) A small magnetic field lifts the degeneracy between states with the same 
value of mp^- A permanent electron electric dipole moment further breaks this degeneracy, but with opposite sign for the 
upper (solid arrow) and lower (dotted arrow) transition. Energy splittings not to scale. 



chopping the sign of the magnetic field amounts to cleanly reversing the direction of time, and will leave certain 
transition energies rigorously unchanged if c?e = 0. These are our "science transitions", which we will measure with 
our highest precision. 

E. Electronic Levels, Spin Preparation, and Spin Readout 

The density of trapped molecular ions will be too low to permit direct detection of the radio frequency or microwave 
science transitions. (A possible exception could involve the use of a superconducting microwave cavity, but this would 
add considerable experimental complexity.) We will of necessity rely on electronic transitions to prepare the initial 
electron spin state, and on a double resonance method to detect the spin flips. The details of these steps will depend 
on the specific molecule we use. For a qualitative illustration, we present a schematic of the calculated l ow- ly ing 
electronic potential curves of HfF"*" (Fig. [5]). We note that HfH+ and ThF+ have similar level structures [sol Il0ll |. 

The molecules will be formed by laser ablation and cooled by supersonic expansion such that a large portion of the 
molecular population will be in ^Sq ground state with a few rotational levels occupied (Sec. Spin-orbit mixing 

between states of identical are enhanced by relativistic effects in the high-Z Hf atom. The b(l) and c(l) states 
are well-mixed combinations of ^ni,'^ni, and ^'Si states, allowing for electric dipole transitions to and from these 
states that do not respect spin selection rules. The ^Sq state, on the other hand, has no nearby = state with 
which to mix, and thus S and A are good quantum numbers. Similarly, the '^Ai state has so little contamination of 
"'^IIi in it th at a rough calculation indicates that it is mctastable against spontaneous decay, with a lifetime of order 
300 ms [30l.fToH. 

The Ramsey resonance experiment will begin with a two-photon, stimulated Raman pulse, off-resonant from the 
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FIG. 2: Potential energy curves for select states of HfF+ [30l ]. The b(l) and c(l) states are well-mixed combinations of ^IIi,"^ Hi, 
and "^E^ states. 



intermediate ^''^IIi states, which will coherently transfer population from the ^Eq, J — Q ground state to the two 
Impl — 3/2 magnetic sublevels of the '^Ai, J = 1 level. The relative phase between the two magnetic levels evolves 
at a rate given by the energy difference. After a variable dwell time, a second Raman pulse is applied, which will 
coherently transfer a fraction of the population back down to the ^Sq state, with probability determined by the 
accumulated relative phase. By varying the dwell time between Raman pulses, the population in the ^Eq state will 
oscillate at a frequency given by the energy difference between the two spin states in the '^Ai manifold. 

The final step in the resonance experiment is to measure the number of molecules remaining in the Ai state. This 
wc propose to do with state-selective photodissociation. Molecules in the ^Ai state will be dissociated via a two-color 
pulse, back up through the "^IIi state to a repulsive curve, generating a Hf+ atomic ion and a neutral fluorine atom. 
Molecules in the ^Eq state will not be affected by the two-color laser pulse and will remain as HfF+ molecular ions. 
The Paul trap parameters will be adjusted to conflne only ions with the Hf"*" atomic mass, and not the HfF+ molecular 
mass with mass difference AM = 19 amu. Finally, the potential on an endcap electrode will be lowered, and the 
remaining ions in the trap will be dumped onto a ion-counting device. 

Details of this procedure will depend on the molecule ultimately selected for this experiment. We are also inves- 
tigating alternative modes of sp in st ate readout, including large-solid-angle collection of laser-induced fluorescence, 
and high finesse optical cavities 

III. EXPERIMENTAL APPARATUS 
A. Molecular Beamline 

We are interested in studying molecular radicals and therefore must create the molecules in situ. As described in 
Sec. Ill Bl we have a small collection of molecules that satisfy our selection criteria and our final choice of molecule 
has not been made. However, for clarity this section will describe the production, detection, and characterization of 
a beam containing neutral HfF molecules and HfF+ molecular ions. 

The molecules are made in a pulsed supersonic expansion (Fig. (H]). A pulse valve isolates ~ 7 atmospheres of argon 
that is seeded with 1% sulfur hexafluoride (SFg) gas from the vacuum chamber. The pulse valve opens for ~ 200 /iS 
allowing the Ar -I- 1% SFg mixture to expand into the vacuum chamber. This creates a gas pulse moving at 550 m/s in 
the laboratory frame, but in the co-moving frame the expansion cools the translational temperature of the Ar atoms 
to a few Kelvin. 

Immediately after entering the vacuum chamber, the gas pulse passes over a Hf metal surface. Neutral Hf atoms 
and Hf+ ions are ablated from this surface with a 50 mJ, 10 ns, 1064 nm Nd:YAG laser pulse. The ablation plume is 
entrained in the Ar -f- 1% SFg gas pulse and the following exothermic chemical reactions occur: 

Hf+SFg HfF + SFg, (1) 

Hf+ + SF6 HfF+ + SF5, (2) 

In the co-moving frame, the resulting neutral HfF molecules and HfF"*" molecular ions are cooled through collisions 
with the Ar gas to rotational, vibrational, and translational temperatures of order a few Kelvin. The molecular 
beam then passes through a skimmer, flrst entering a region where laser induced fluorescence (LIF) spectroscopy is 
performed and finally arriving at an rf (Paul) trap where the ions are stopped and confined. 
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FIG. 3: Experimental setup. Laser ablation of a metal Hf target creates neutral Hf atoms and Hf*" ions that react with SFe to 
produce neutral HfF molecules and HfF^ molecular ions, respectively (Eqs. [1] and[2]). The molecules (both neutral and ionic) 
are cooled in a supersonic expansion with a He bufTer gas. The molecular beam is illuminated with a pulse dye laser beam 
and the resulting fluorescence is collected with a photomultiplier tube (PMT) yielding laser induced fluorescence (LIP) spectra 
(Fig. |4]). At the end of the beamline, the ions are loaded into an rf (Paul) trap where the electron spin resonance experiment 
is performed. The Paul trap also acts as a quadrupole mass filter and ions of a particular mass/charge ratio are detected with 
a microchannel plate (MCP) (Fig. [5]). Additionally, the spatial resolution of the MCP allows for the temperature of the ion 
cloud to be determined from the detected cloud size. 




14216 14220 14224 14228 

photon wavenumber [cm"^] 



FIG. 4: Laser induced fluorescence (LIF) spectroscopy. The top trace is experimental data for a newly detected neutral HfF 
transition: [14.2] \Q\ — 3/2 \v' = v" , J') <— X^A3/2 \i^",J")- The transition highlighted with a vertical arrow originates from 
the rotational ground state. The bottom trace is a theoretical prediction assuming a rotational temperature of 5 K. The traces 
are offset vertically for clarity. 



LIF spectroscopy is performed by transversely illuminating the molecular beam with a 500 /iJ, 10 ns, ~ 700 nm 
dye laser pulse. The linewidth of the dye laser is specified to be less than 0.1 cm~^. Fluorescence photons are collected 
and imaged onto a photomultiplier tube (PMT). 

Using this technique we have found previously unobser ved neutral HfF molecular transitions, one of which is shown 
in Fig. m (for previous neutral HfF spectroscopy see Ref. |ll2| ). The data shows that entrained neutral HfF molecules 
are cooled to rotational temperatures of order 5 K, with a large fraction of the population in the rotational ground 
state. We expect that entrained HfF+ molecular ions should be similarly cooled. 

To detect the presence of HfF+ molecular ions in the beam the rf (Paul) trap is operated as a quadrupole mass 
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FIG. 5: Mass spectrometry. Operating the rf (Paul) trap as a quadrupole mass filter gives mass-dependent trapping potentials 
such that Hf+ (M = 180 amu), HfF+ (M = 199 amu), H£F+ (M = 218 amu), and HfF+ (M = 237 amu) can be separately 
trapped and detected. The ion detector signal is a non-linear function of ion number, but a level of 0.4 corresponds to ~ 100, 000 
ions. 



filter. All of the ions in the beam are stopped and loaded into the trap. The voltages applied to the trap electrodes 
are then adjusted only to confine ions of a particular mass/charge ratio. Finally, the ions remaining in the trap are 
released onto the ion detector and counted. A typical mass spectrum is shown in Fig. El which clearly resolves the 
HfF"*" molecular ions from the other atomic and molecular ions in the trap. 

Our experimental count rate will be limited by space charge effects of the trapped ions. Therefore, any ions trapped 
that are not used in measuring the eEDM limit the statistical sensitivity of our measurement. In order to maximize 
our count rate, we wish to create and trap only HfF'*' ions of a single Hf isotope and in a single internal quantum 
state. One scheme is to filter out all of the ions created from laser ablation and use photoionization techniques to 
ionize neutral HfF in as state-selective a way as possible. Using two color, two photon excitation, we e xcite to a high 
lying Rydberg state, in an excited vibrational level, that then undergoes vibrational autoionization [ll3l |. The ion 
core of these Rydberg state molecules will occupy a single rotational level and consist of a single Hf isotope. The 
autoionization process is seen, in our preliminary (unpublished) data, to leave the ion core rotational level largely 
unperturbed. It should be possible to excite a Rydberg level that corresponds to an excited '^Ai ion core with v = 1, 
J = 1 (where v is the vibrational quantum number). The Rydberg state might then vibrationally autoionize to the 

= 0, J = 1 -^Ai level that will be used to measure the eEDM. 



B. Radio Frequency (Paul) Trap 

For our preliminary studies of ion production, the ions are confined by a linear rf (Paul) trap shown schematically in 
Fig. El The ideal hyperbolic electrodes are replaced by cylinders of radius aw 1.15po, where po is the minimum radial 
separation between the trap center and the surface of the electrodes. This choice produces the best approximation to 
a perfect radial two-dimensional electric quadrupole field [115]. 

For HfF"*" (M = 199 amu), an example set of operating parameters for the ion trap would be po — 25 mm, 
Vri = 550 mV, and Wrf = 27r x 15 kHz. This produces a ponderomotive potential that is well within the harmonic 
pseudo-potential approximation given by J7rf(p) = Moj'^^^p^ /2, where the radial secular frequency is approximately 
i^scc = with q — AeVy^f/Mpluj^f. For the above parameters, q = 0.2, Wsoc = 2n x 1 kHz, and J7rf(po) — 300 K. 

Under these conditions, an ion cl oud at a temperature of 15 K would have an rnis radius of 5 mm. The trap can also 
be operated in mass filter mode [ll6j |. 

In addition to supplying the oscillating electric quadrupole field for radial confinement, the cylindrical electrodes 
can also be driven with voltages to produce the rotating electric bias field, frot, needed to polarize the molecular ions 
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FIG. 6: Linear rf (Paul) trap. Neighboring cylindrical electrodes are driven with rf voltages 180° out of phase. Axial confinement 
is provided by d.c. voltages applied to the end cap electrodes. The cylin drica l electrode rods have radius a and the radial distance 
from the trap center to the nearest electrode surface is po- See Ref. [ll4( | for further details of rf (Paul) trap operation. In 
addition to the voltages oscillating at LUrf, there is also a component of the voltages oscillating at Urot. Over a period of time 
27r/a;rot, the electric field at the a^dal center (z=0) of the trap will trace out a trajectory which subtends a solid angle A of 
exactly 2n. Ions to the left (right) of trap center will experience an electric field whose trajectory subtends slightly less (greater 
than) 2n. Consequences of this time variation are explored in discussed in Sec. IIV Cl and HV Dl Not to scale. 



(Fig. in]). In order to generate £iot neighboring electrodes will be driven 90° out of phase at a frequency Wrot- The net 
voltage applied to each electrode is the sum of the voltages Vif + V^ot • 

At present we are designing a second-generation ion trap with geometry designed for optimal precision eEDM 
spectroscopy, rather than for mass selection. The perfect ion trap would have very large optical access for collection 
of laser-induced fluorescence, and idealized electric and magnetic fields as follows 

£ = £rotP + £^fixx - yy) cos(wrfO 4- £'A-zz + yy/2 + xx/2) (3) 



B = Brotp' (4) 

where p' = cos(a;rotO^ + sin(a;rotOy ^n'i ^rf ~ ~2Vrf/po- 

If we assume Wrot ^ Wrf, that Wrot/wrf is not a rational fraction, and that w^f ^ e£^f/M, then we can cleanly 
separate out the ion motion into three components: rf micromotion, circular micromotion, and secular motion. 

rf micromotion involves a rapid oscillation at u-^ whose amplitude grows as the ion's secular trajectory takes it 
away from trap center. The kinetic energy of this motion, averaged over an rf cycle, is given by 



where x and y in this case refer to the displacement of the ion's secular motion. 
The displacement of the ion's circular micromotion is given by 



—■ (6) 



The kinetic energy of the circular motion, averaged over a rotation cycle, is given by 
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The time-averaged kinetic energies of the two micromotions act as ponderomotive potentials that contribute to the 
potential that determines the relatively slowly varying secular motion: 

Usee = E,ot (2;, y, z) + Erf(x, y, z) + eS',{2z^ - if - x^)/A. (8) 

In the idealized case, the secular motion corresponds to 3-d harmonic confinement with secular or "confining" fre- 
quencies 



1 d^Us, ^ 



^SCC 



M di^ 



(9) 



for i = x,y, z. In the idealized case, confinement is cylindrically symmetric, uJx = ^y, and S^ot is spatially uniform, so 
the circular micromotion does not contribute to the confining frequencies. 

The density of ions will be low enough that there will be few momentum-changing collisions during a single mea- 
surement. Thus, any given ion's trajectory will be well approximated by the simple sum of three contributions: 

(i) a 3-d sinusoidal secular motion, specified by a magnitude and initial phase for each of the x, y, and z directions. 
In a thermal ensemble of ions, the distribution of initial phases will be random and the magnitudes, Maxwell- 
Boltzmannian. For typical experimental parameters (see the Appendix) the secular frequencies Wi will each be about 
27r X 1 kHz and the typical magnitude of motions, r, will be about 0.5 cm. 

(ii) the more rapid, smaller amplitude rf micromotion, of characteristic frequency about 2tt x 15 kHz and radius per- 
haps 0.05 cm. This rf micromotion, purely in the x-y plane, is strongly modulated by the instantaneous displacement 
of the secular motion in the x-y plane, and vanishes at secular displacement x=y=0. 

(iii) The still more rapid rotational micromotion, purely circular motion in the x-y plane, at frequency Wiot about 
27r X 100 kHz and of radius comparable to the rf motion, around 0.05 cm. In the idealized case, the rotational 
micromotion (in contrast to the rf micromotion) is not modulated by the secular motion. 

As described in Sees. ITVEl and Ivl below, for spectroscopic reasons we must operate with trapping parameters such 
that Ej-ot ^ SOksT. Under that condition, relatively small imperfections in £rotj say a spatial variation of 1.5%, can 
give rise to contributions to Usee of the same scale as the ions' thermal energy, and thus significantly distort the shape 
of the trapped ion cloud or even deconfine the ions. 

For improved optical access we had to shrink the radius of the linear electrodes a with respect to their spacing po 
The spectroscopic requirement for highly uniform fiot then forced the redesign of the second-generation ion trap to 
be based on six near-linear elements arranged on a hexagon, rather the four electrodes arranged on square shown in 
Fig. [HI The trap will be discussed in more detail in a future publication, but simulations project spatial uniformity of 
£iot better than 0.5% with good optical access. The design led to significant compromises in the spatial uniformity of 
£rf, so in future operation, mass selectivity in ion detection will come not from a quadrupole mass filter, but rather 
from pulsing frot to a very high value for a small fraction of a rotation cycle and then doing time-of-flight mass 
discrimination on the ions thus ejected. B^ot will be imposed by means of time-varying currents flowing lengthwise 
along the same electrodes that generate frot- 



IV. SPECTROSCOPY IN ROTATING AND TRAPPING FIELDS 



On the face of it, an ion trap, with its inhomogeneous and rapidly time-varying electric fields, is not necessarily 
a promising environment in which to perform sub-Hertz spectroscopic measurements on a polar molecule. In this 
section we will explore in more detail the effects of the various components of the electric and magnetic fields on 
the transition energies relevant to our science goals. The theoretical determination of the energy levels of heavy 
diatomic molecules in the presence of time-varying electric and magnetic fields is a tremendously involved problem 
in relativistic few-body quantum mechanics. State-of-the-art ab initio molecular structure calculations are limited to 
an energy accuracy of perhaps 10^'^ Hz, a quantity which could be compared with the size of a hypothetical "science 
signal" , which could be on the order of 10^'' Hz or smaller. 

Fortunately, we can take advantage of the fact that at the energy scales of molecular physics, time-reversal invariance 
is an exact symmetry except to the extent that there is a time-violating moment associated with the electron (or 
nuclear) spin. In this section, except in those terms explicitly involving d^, we will assume that time-reversal invariance 
is a perfect symmetry in order to analyze how various laboratory effects can cause decoherence or systematic shifts 
in the relevant resonance measurements. The results can be compared to the size of the line shift that would arise 
from a given value of the electron EDM, which is treated theoretically as a very small first-order perturbation on the 
otherwise T-symmetric system. 

In the subsections below, we bring in sequentially more realistic features of the trapping fields. 
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A. Basic Molecular Structure 



We begin by considering in detail the relevant molecular structure in zero electric and magnetic fields, thus quan- 
tifying the qualitative discussion of the experiment given in Sec. [TTl Although the molecular structure cannot be 
calculated in detail from ab initio structure calculations, nevertheless its analytic structure is well known. Because 
the measurements will take place in nominally a single electronic, vibrational, and rotational state, we will employ an 
effective Hamiltonian within this state, as elaborated by Brown and Carrington jll7l |. This approach will specify a 
few undetermined numerical coefficients, whose values can be approximated from perturbation theory, but which will 
ultimatel y be meas ured . 

Brown [TT8| - [T20| and co-workers have done thorough work on deriving an effective Hamiltonian for molecules. 
The complete Hamiltonian in the absence of de is given by 

instruct = -ffcloc + -ffvib + Hso + Htum + i?SS + -^SR + i?HFS + i?LD , (10) 

listed in rough order of decreasing magnitude. Since we are concerned only with terms acting within the subspace of 
the '^A manifold, other electronic and vibrational states will enter only as perturbations that help to determine the 
effective Hamiltonian. Thus we consider eigenstates of iJeiec and i?vib- 

The remaining terms in Eq. (jlOl) are corrections to the Born-Oppenheimer curves. They describe couplings between 
various angular momenta (i?sRi H^ifs), parity splittings {Hi^u, H^ifs), and spin-dipolar interactions (i?ss, -f^HPs)- In 
typical Hund's case (a) molecules these interactions are small compared to the rotational energy governed by i?tum- 
The relevant interactions that act within the = 1 manifold of states take the explicit form 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

The constants in the first four terms are as follows: A is the molecular spin-orbit constant. Be the rotational constant 
for the electronic level of interest, D the effect of centrifugal distortion on rotation (typically D ^ Be{me/myao\)^ , 
with rrie the electron mass and mmoi the reduced mass of the molecule), A governs the strength of the spin-spin 
dipolar interaction, and 7sr determines the strength of the interaction of the spin with the end-over-end rotation of 
the molecule. These four terms primarily describe an overall shift of the '^Ai J-level, and can be ignored in evaluating 
energy differences in the states we care about. They can, however, contribute small perturbations to these basic levels, 
as we will describe below. 

Within the '^Ai, J = 1 manifold of interest, the energy levels are distinguished by the hyperfine and A-doubling 
terms. The hyperfine Hamiltonian Hufs includes the familiar contact (bp), nuclear-spin-orbit (a) and spin- nuclear 
spin terms (c). By estimating the parameters in perturbation theory, it is expected that the resulting hyperfine 
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splitting is on the order of 27r x 50 MHz 10l| . The hyperfine interaction also contains a previously unreported term, 
with constant denoted ca, that is connected to the A-doubling. This term is expected to be even smaller than the 
already small A-doublet splitting itself |l2l| . however, and will be ignored. 

The A-doubling Hamiltonian arises from Coriolis-type mixing of states with differing signs of A due to end-over-end 
rotation of the molecule. For a "^A state this interaction is characterized by three constants, of which the parameter 
OA is the dominant one. These terms describe how the ^A state is perturbed by electronic states with ^S-i-ijj ^^^^ 
^^"""^E symmetry. Since we are primarily concerned with terms in the Hamiltonian that affect the ground rotational 
state of the "^Ai electronic level, we only need to keep the term which connects r2 = ltoSl = --l. This term has 
the general form, with numerical prefactors Cn,s,n' that depend on Clebsch-Gordon coefficients and wavefunction 
overlap, |118i] 

|0A + 3pA + 6,A| = 5a ^ ^n^-^-(EA-En)(EA-''E;)(EA-EnO' ^''^ 



where the sum is over all intermediate E and H states of singlet and triplet spin symmetries. For HfF"*" this perturbation 
leads to a A-doublet splitting on the order of 27r x 10 kHz. This estimate was carried out assuming a aS molecular 
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orbital configuration, where the 5 orbital has total angular momentum L = 2 in the pure precession approximation. 
The ground X^H is a molecular orbital but has some admixture of atomic do orbitals. We therefore expand the 
molecular wavefunction into atomic orbitals and reduce the amount of admixture by the factor Cd that describes the 
do character. From here on, we shall express the energy difference in parity levels for the J = 1 as iOd — 46a, rather 
than 6a itself. 

Thus the basic molecular structure of interest to the '^Ai, J = 1 state is governed by two constants: the hyperfine 
splitting Ehf (given by 3A||/4 for J = 1,7 = 1/2) and the A-doublet splitting Wef- These constants give the structure 
depicted in Fig. [Ua). These basic levels may be perturbed by couplings to other levels, especially rotational or 
electronic excited states. However, for the J = 1 state of interest, some simplifications are possible, namely: (1) 
Off-diagonal couplings in 51 are zero since J • S preserves the value of J (there is no level with J = 1 and 51 = 2); 
(2) Off-diagonal contributions that mix J = 2 into the J — \ manifold thus depend solely on the applied fields and 
the hyperfine interactions. Since the value of the spin-orbit constant is expected to be far larger than the rotational 
constant and we are concerned with a J = 1 state, the operators that connect 51 to 51 ± 1 will be ignored. The 
contributions to the ground state characteristics by terms off diagonal in 51 are smaller by a factor of the hyperfine 
interaction energy to the spin-orbit separation energy, hence a factor of 10^^. This is the value which appears in front 
of any term connecting 51 to 51 ± 1 in the ground J = 1 state. 



B. Effect of Non-rotating Electric and Magnetic Fields 

The influence of external fields presents new terms in the Hamiltonian of the form 

-ffstark = — dmf • £ (18) 
T^Zccman = -[l-B. (19) 

Here £ and B are the electric and magnetic fields, assumed for the moment to be collinear so that they define the 
axis along which mi? is a good quantum number; while dmf and \i are the electric and magnetic dipole moments of 
the molecule. 

The electric dipole moment arises from the body- fixed molecular dipole moment, at fields sufficiently small not 
to disturb the electronic structure. We assume that the field is sufficiently large to completely polarize this dipole 
moment, i.e., dmff ^ Wef, in which case the Stark energies are given by 

Estark = -mpnjpdn^f£, (20) 

where 7f is a geometric factor, analogous to a Lande g-factor, which accounts for the Stark effect in the total angular 
momentum basis F. In the limit where the electric field is weak compared to rotational splittings, it is given by 

_ J(J + l) + f(F+l)-/(/ + l) 

2^^(F + 1)J(J+1) ■ ^ ' 

Its numerical values in the J = 1 state are therefore 7_f=3/2 — 1/3 and 'yF=i/2 = 2/3. The electric field therefore 
raises the energy of the states with mpn < (denoted "upper" states with superscript u), and lowers the energy of 
states with itif^ > ("lower" states with superscript £). This shift in energy levels is shown in Fig. [ijb), where \a) 
and \b) are upper and |c) and \d) are lower states. 

The form of the Zeeman interaction is somewhat more elaborate, as the magnetic moment of the molecule can arise 
from any of the angular momenta L, S, J, and I. Quite generally, however, in the weak-field limit where (j-bB <^ Ehf, 
the Zeeman energies are given by rriFgp^^ ^bB, where ^b is the Bohr magneton and gp^^ are g-factors for the upper 
and lower states. In general, gp ^ gp, and this difference can depend on electric field, a possible source of systematic 
error. We will discuss this in Sec. lIVGl below. 

The leading order terms in the Zeeman energy are those that preserve the signed value of 51. They are given by 

-ffzoGinan = {jF [((.9L + 9r)^ + (QS + .9r)S) 51 - J (J + 1)] - giKp) Ulp^iBB, (22) 

where up = {F{F + 1) + /(/ + 1) — J(J + 1))/2F{F + 1) is another Lande-type g-factor, but for nuclear spin. The 
orbital and spin g-factors are gp and gs, while the rotation and nuclear spin g-factors are gr and gj. Both gr and gi 
are small, being on the order of the electron-to-molecular mass ratio ~ nie/m^oi ^ 10~^. Thus for an idealized '^Ai 
molecule where gp — 1, A — ±2, gs — 2, S — =f1, we would expect molecular g-factors on the order of 10~^. More 
realistically, gs differs from 2 by a number on the order of a, the fine structure constant, and a g-factor ^ 10~^ might 
be expected. In heavy-atom molecules such as ours for which spin-orbit effects mix A, we may expect instead the 
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difference 2gL — gs to be as large as ~ 0.1 in magnitude. If we assume the dominant contribution comes from these 
spin-orbit type effects, we can define the (j-factor for the J = 1 state as 

5F=3/2 = 7f=3/2(.9lA + gs^)n < 0.03, (23) 

while 

9f=i/2 = 25f=3/2- (24) 
Finally, the effect of the EDM itself introduces a small energy shift 

HeDM = -de ■ £cS — de£cS<^i • n, (25) 

where ai is the spin of the s-electron contributing to the EDM signal; and fi denotes the intermolecular axis, with n 
pointing from the more negative atom to the more positive one; in our case from the fluorine or hydrogen to thorium, 
platinum, or hafnium. Also in this convention we take fcff as positive if it is anti-parallel to h. The energy shift 
arising from this Hamiltonian depends only on the relative direction of the electron spin and the internuclear axis, 
and is given by 

Polarizing the molecule in the external field selects a definite value of fi, hence a definite energy shift, positive or 
negative, due to the EDM. This additional shift is illustrated in Fig. [TJc). 

For a range of field strengths and parameters, the energies of the sublevels within the J = 1 manifold are well 
approximated by a linear expansion in the electric and magnetic fields. We define 

S = S||^+6_L (27) 

Taking ^ di^^iE <C Ehf and dmf^ 3> .gFMS'8||, and setting B±^ = 0, we get for the non-rotating energies, 

W^l'{F,mF,^;£,B) « ^{F{F + l) - ^)^M-mpn-iFd,,,i£ + mpgl''^iBB-{d,£,fi/2M)n, (28) 

where Q, is either 1 or -1, and the prefactor in front of Ehf is such that for the J — 1 level, E(i^ = 3/2) - E(F — 1/2) 
= 3A||/4 = Ehf. F and £7 are good quantum numbers only to the extent that the electric field is neither too large 
nor too small, but we will use F and as labels for levels even as these approximations begin to break down. 
For notational compactness, we introduce special labels for particular states as follows (see Fig. (Hb)): 

(29) 



\a) = 




3/2, m = 


3/2, = -1) 


\b) = 




3/2, m = 


-3/2, f7 = 1) 


\c) = 




3/2, m = 


3/2, = 1) 


\d) - 


\F = 


3/2, m = 


-3/2, f7 = -1 



with corresponding energies, Eq, Ef,, Ec, and E(i, and identify the energies of two particularly interesting transitions, 
W^" = Ea - Eh, and = E^ - E^ such that 

= -igl^iBB + d,e,s 
W' = ig'p^iBB - d,£,s. (30) 

Taking this analysis a step farther, it is possible that the electric field energy dmfS is not small compared to the 
hyperfinc splitting Ehf. In this case the electric field mixes the different total- states and perturbs the above energies. 
Ignoring the magnetic field and EDM energies, the energy levels take the form 



E„,.(F - 3/2, mp^ = +1/2) 



1 / d^fS 
2 



E„,(F^ l/2,mFf7 = +l/2) = -- 



1 ( dynfS 

2 



Ehf 
~3~ 

Ehf 

~3~ 



dmt£ 



Ehf 



d-^f£ 



dm{£ 



E 



hf 



dnifS 



3^gf 

4dmff 
Sdniff 



(31) 
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EnriF ~ 3/2, m^f^ = -1/2) 



Enr(F ~ 1/2, mF^i = -1/2) 

The equations of this section have so far been to one degree or another approximate results. But in the absence of 
exotic particle physics we can invoke time-reversal symmetry and write exact relations: 

E„,(F, mp, n;£,B) - E„r(F, -m^, -f^; f , B) = E„,(i^, -mp^-n; f , -B) - Enr(F, mp^n- £, -B) (33) 

which, for B — 0, becomes 

E„,(F,mF,f^;f) = E„,(F,-mF,-rj;f). (34) 

This exac t degeneracy is, in fact, an example of the Kramers degeneracy that follows from time-reversal invari- 
ance |l22l |. For our purposes, the key result here is that, in the limit of non-rotating fields, zero applied magnetic field, 
and an electron EDM, the energy of the science transitions \mF,^) ^ | — mp,—^) (and in particular, and W^) 
are independent of the magnitude of the electric field. This is an important property because we are using spatially 
inhomogeneous electric fields to confine the ions in the trap, and we want to minimize the resulting decoherence. 

This degeneracy in turn means that the energy differences and W'' depend only on the magnetic field and, of 
course, the EDM term as shown in Eq. pOI) . The magnetic contribution reverses sign upon reversing the direction of 
B with respect to the electric field direction (which also sets the quantization axis, since dmfS ^ fJ^sB). Therefore 
the science measurement is given by the combinations 

Wi£,B) + W{£,-B) = 2d,£,s 

W\£,B) + W\£,-B) = -2de5cff, (35) 
where a -I- sign on B denotes that it points in the same direction as £ . 



1 / di^i£ Eht 
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dmf£ 
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2 o 2 
cf 



(32) 



ai£ ' 



C. Rotating Fields, Small-Angle Limit 

Many EDM experiments over the years have been complicated by the problem of "Berry's phase" , the term in this 
context used as a catch-all to describe a variety of effects related to the motion of the particles in inhomogeneous 
fields. 

The sketch in Fig. [Tlja) illustrates the classic Berry's phase result: if the field that defines the quantization axis, as 
experienced locally by a particle (or atom, or molecule), precesses about the laboratory axis at some angle, 9, then, 
in the limit of slow precession, with each cycle of the precession the wave- function 5* picks up a phase given by mpA, 
where is the instantaneous projection of the particle's total angular momentum on the quantization axis, and A 
is the solid angle subtended by the cone. If the precession is periodic with period r, one can (with provisos, as we will 
discuss) think of this phase-shift as being associated with a frequency, or indeed energy, mpA/r . In a spectroscopic 
measurement of the energy difference between two states whose mp values differ by Smp, there will be a contribution 
to the transition angular frequency ASmp/r. 

In neutron EDM experiments, motional magnetic fields, in combination with uncharacterized fixed gradients from 
magnetic im purit ies, Berry's phase can be a dangerous systematic whose dependence on applied fields can mimic an 
EDM signal jl23j |. In Sec. lIVLl we will see that the effects of motional fields in our experiment are negligible. 

Neutral atoms or molecules may be confined in traps consisting of static configurations of electric or magnetic 
fields. These traps are based on the interaction between the trapped species' magnetic or electric dipoles and the 
inhomogeneous magnetic or electric fields, respectively, of the trap. Especially in cases where the traps are axially 
symmetric, so that the single-particle trajectory of an atom can orbit many times one way or the other about the axis 
of the tra p, th e coherence time of an ensemble of atoms with a thermal distribution of trajectories can be severely 
restricted [l2J|. Our system is quite different, because in an ion trap the forces arise from the interaction between 
the trapping fields and the monopole moment of our trapped ion. Assuming the temperature, size of bias field, and 
radius of confinement are the same, the trapping fields for an ion are spatially much more homogeneous than would 
be those for a neutral molecule or atom. 

That said, the fact that we can speak of a "bias" electric field at all in an ion trap conies at the cost of having the 
applied electric field constantly rotating. 
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(b) 



FIG. 7: (a) Small-angle limit. When the quantization axis T follows a slow periodic perturbation characterized by tilt angle 
6^, angular frequency oj and enclosed solid angle A, two states whose instantaneous projection of angular momentum along T 
differs by 5m will have their effective relative energy displaced by a "Berry's energ" loAimli-K. (b) Large-angle limit. When 
instead the quantization axis sweeps out a full 2-k steradians per cycle (a^O), the differential phase shift between the two levels 
is indistinguishable from zero, and in the most natural conceptual framework, the Berry's energy vanishes. 



D. Rotating Fields, Large-Angle Limit (Dressed States) 

The basic dressed-state idea is an extension of the more common idea of an energy eigenstate: a system governed 
by a time-invariant Hamiltonian H will have solutions 5* such that \E'(i -|- T) = e~'"^\E'(t) for all T and t\ such a 
solution 5* is called an energy eigenstate, with uj being then the corresponding energy. Similarly, a system governed by 
a periodic Hamiltonian with period t such that H{t + T) = H{t) for all values of t, will have so-called "dressed-state" 
solutions 4* such that 'ii{t + nr) — e"™'^4'(i) for all t and all integer values of n. It is tempting to call 0/t the 
"energy" of the dressed state, but there will be an ambiguity in that energy because we can always replace </> with 

Operationally, the dressed state energies are derived from the eigenvalues of a formally time- independent Hamilto- 
nian. If Hq denotes the Hamiltonian in the absence of the field, then the appropriate rotation-dressed Hamiltonian is 
given by 

J^drcsscd = Hf) — dmf ' i^^rot + ^^rot, (36) 

Hiot is defined as [125 1 

Hrot = -uJrot {cos{9)F, - sm{e)F,) (37) 

where i^^^and Fj. are the projections of the total angular momentum F into a set of axes where z coincides with the 
instantaneous direction of the electric field. We now make explicit the rotating electric field with frot- The cos(6') 
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term thus provides an energy which, when multiphed by the rotational period r = 27r/a;iot, gives the ordinary Berry 
phase, 

- 27rcos(6')mF ^ 27r(l - cos(6'))mF (38) 

where we have taken the hberty of adding an arbitrary phase 2'KmF to reveal explicitly the solid angle 27r(l — cosfl). 

In the experiment, the applied electric field should lie very nearly in the plane orthogonal to the rotation axis, i.e., 
9 K, 7r/2. It is therefore useful to consider the small angular deviation from this plane, a = -k 12 — (Fig. [7]). Then 
the apparent energy shift arising from the geometric phase is 

Egco = -mpuj^ot su\{a) « -mpu^ata. (39) 

Now consider two states which are, in the absence of rotation, degenerate, say the states |a), with m = 3/2, f7 = — 1, 
and the state with m — —3/2,17 = 1, indicated in Fig[ljb). Rotation breaks this degeneracy, by adding the 
energies ~ ±[2>/2)u)rota^ as shown by the dashed hues in Fig.[8l These levels cross at a = 0, leading to their apparent 
degeneracy when the electric field lies in the horizontal plane. 

In addition, the rotation of the field also incurs coupling between states with different mp values, arising from the 
sin(0) term in Eq. (j37[) . This perturbation, treated at third-order in perturbation theory, connects the two levels and 
turns the crossing into an avoided one, as shown by the solid lines in Fig. [51 Since the energy contribution due to the 
rotating field is small compared to Stark energy splittings, we can use ideas similar to the derivation of A-doubling, 
i.e., we take a sum of the perturbing components and take them to the appropriate power. We look for terms in this 
expansion that can connect the state \a) = \mp^) to 1^*) = | — mp ~ fi). Therefore, the power of perturbation theory 
needed is 2mF + 1, where the 2mF takes mp — > —mp and the extra power takes O — — f2. The two terms in the 
Hamiltonian that can do this are the A-doubling term and the mj^^-changing terms of the rotating electric field. Our 
expansion is, schematically, the following 



(g LD+grot)''"^ + ^ 
(AE„,J2mp 



^coup- Tx-E^ — ^2;;;;^ • (40) 



The (AEnip)^™'' are the energy level differences between states with different mp values, thus are related to the Stark 
splittings. This tells us that 

/ \ 2mF 

A-^.(-^) , (41) 

where A is the energy splitting at the level crossing between otherwise degenerate states with mp > Q and mp < Q. 
The numerical prefactor in this expression has a rather complicated form within perturbation theory. Howev er, it s 
value can be computed by numerically diagonalizing the relevant hyperfine-plus-rotation dressed Hamiltonian (121| . 
The result, for the mp = ±3/2 states in Fig.[TJb), is 

^-/^^nQ^,J-^\\ (42) 

where the superscript u refers to mixing between the \a) and \b) states, and the superscript t to mixing between |c) 
and \d) states. In the absence of the hyperfine interaction, the average value of the numerical prefactor is 170 and 
the upper and lower states have the same avoided crossing. However, small fractional differences between A" and A^ 
turn out to be significant, and are discussed further below. 

The presence of the electric field causes the states with \F = 1/2, mp = ±1/2) and \F — 3/2, m^? = ±1/2) to mix. 
Including the hyperfine interaction into the numerical diagonalization yields 

A = \{^- + AO - 170a.ef (:r^)\ (43) 



<5A=-(A"-A0 « 127..ef Mp^ • (44) 

It is evident that the average shift is the same, but now the upper and lower levels acquire a different splitting due 
to the rotation-induced mixing within the sublevels. The difference is suppressed relative to the average value of the 
splitting by a factor of {d„^f£ /Ei^f)'^ , reflecting the fact that higher orders of perturbation theory are needed to include 
the effects of the hyperfine interaction. For oj^f = 27r x 10 kHz, Wj-ot = 2% x 100 kHz, o?mf£i.ot = 27r x 10 MHz, E^f = 
27r X 45 MHz, then A = 27r x 2 Hz and 5A = 27r x 0.06 Hz. 
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FIG. 8: The apparent energy shifts between mp — +3/2 and mp = —3/2 states in upper (a,b) and lower (c,d) A-doublet levels 
versus a, the angle of the electric field to the plane orthogonal the rotation axis of Srot (a is shown in Fig. [Zib)). (a) At a = 0, 
there is an avoided crossing that mixes mp = ±3/2 states, with an energy splitting at the crossing of A"/^ (b) Since a = at 
the axial trap center, and since we need ttif to be a signed quantity in order to measure de, we will bias away from the avoided 
crossing using a magnetic field Srot. SmpOF f^sBrot > A"'^* is required for mp to be a quantity of definite sign. This picture 
is intuitively correct in the limit that A"'* > ojmax (see Sec. IIVE|| . The experiment will be performed in the opposite limit. 
However, solving the time dependent Schrodinger equation (Eq. [55} gives the same requirement of SmpgE ^J^BB^:ot > A"/^ in 
both limits. 



The magnitude of the rotation- induced mixing within any of the four pairs of otherwise degenerate m — ±1/2 states 
is much larger than the mixing within either pair of m = ±3/2 states, A" or A^. For this reason, the m — ±1/2 levels 
are probably not great candidates for precision metrology in rotating fields. 

An ion in a trap will feel an axial force pushing it towards the axial position where the axial electric field vanishes, 
that is, the location at which a is identically zero. This poses a problem, because at a = 0, each dressed state is an 
equal mixture of states with = 1 and with f2 = — 1. In other words, the dressed states right at the avoided crossing 
will have vanishing eEDM signal. The solution is to bias the avoided crossing away from a = by adding to the 
trapping fields a uniform, rotating magnetic field which is instantaneously always parallel or anti-parallel to vector 

^rot ! 

B,ot = ^£rot- (45) 

In our convention, S^ot defines the quantization axis, so that the number fj-ot will always be taken to be positive. The 
sign of Slot then determines whether the co-rotating magnetic field is parallel {Bmt > 0) or anti-parallel (Srot < 0). 

The energy levels are now as shown in Fig. [5Jb). As derived below in Sec. lIV El in the limit BiotgFfJ-B ^ A"/^, the 
dressed states near a = are once again states of good mp and fi. The energy splitting between the two states, as 
altered by the rotation of the field, are given approximately by 

/e (A"/^)^ 
where the + sign corresponds to u states, and the — sign to i states. 
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Over the course of one axial oscillation of the ion in the trap, a which is approximately proportional to the axial 
electric field, will average to zero. Unfortunately, the contributions to 5W from E^ot and from B^ot are larger than 
that from the scale of the physics we most care about, de^eff, and the spatial and temporal variation in E-^ot and in 
Brot will reduce the coherence time of the spectroscopy, as discussed in Sec. IIV Hi - IIV Jl below. But to the extent that 
one is able quite precisely to chop B^ot to —B-cot on alternate measurements, the science signal still arises from the 
same combination as in Eq. 1351 

{W^"{E.ouB,ot)) + {W^"{E.ou -B,ot)) - ±24fcff , (47) 

where the +/- corresponds to the u/£ superscripts respectively, and the brackets denote averaging over the excursions 
of a, which is assumed to vary symmetrically about zero. 

The equation above relies on several approximations. One needs in particular that dniffj-ot ^ ^cf: S^F/is-Brot ^ 
ITOwof (wj-ot/c^mf^rot)'^, and dmt^rot ^ ^rot, Q! ^ 1, and dmfi^^rot < Ehf- Thcsc are all good approximations, but they 
are not perfect. For example, using values from the Appendix, a;rot/('^infiS'rot) ~ 0.01, a small number, but not zero. 
To what extent will imperfections in these approximations mimic an eEDM signal? 

The driving principle of our experimental design is to measure de with as close to a null background as possible. 
We are not especially concerned if the right hand side of Eq. (jT7|) is 1.9 deScS rather than 2.0 deScS- More important 
to us is that, if de = 0, the right-hand side of Eq. (|T7)) be as close to zero as possible. As we shall see, as long as 
we preserve certain symmetries of the system we are guaranteed a very high quality null. A preliminary remark is 
that the "energy" of a dressed state, or more precisely the phase shift per period r, is unaffected by an offset in how 
the zero of time is defined. A second observation is that, in the absence of exotic particle physics (such as nonzero 
eEDM) , the energy levels of a diatomic molecule in external electromagnetic fields are not affected by a global parity 
inversion. 

Under the action of this inversion, all the fields and interactions in the Hamiltonian transform according to their 
classical prescriptions, whereas quantum states are transformed into their parity-related partners. In a parity-invariant 
system, parity thus changes quantum numbers, but leaves energies of the eigenstates unchanged. This is true for the 
dressed states as well, since their eigen-energies emerge formally from a time-independent Hamiltonian. 

To formulate the effect of inversion symmetry we write the electric and magnetic fields as 

£ = £,otP + £zz (48) 
B = 6„t/5' (49) 

where p' = cos(aj,.oti)^ + sin(a;roti)y and a = tan^^(£2/£'i.ot)- The dressed states defined by the rotating field are 
characterized by the projection of total angular momentum on the axis defined by the rotating electric field, 
^rot/^rot- Because the magnetic field is not strictly collinear with the electric field, and because of the field rotation, 
mp is only approximately a good quantum number. Nevertheless, considering the effect of parity on all the m^'s 
simultaneously, we can still map each dressed eigenstate into its parity-reversed partner. 

Assuming the ions are "nailed down" in their axial oscillation, at a particular value of £z and thus a, our various 
spectroscopic measurements would give dressed energy differences E(£i.ot, Srot, a, m^, ^) — E(£i.ot, ^rot, a, —mp, — fi). 
Now we invoke the following symmetry argument: if we take the entire system, electric fields, magnetic field, and 
molecule, and apply a parity inversion, that will leave the energy of the corresponding levels unchanged. If further we 
then shift the zero of time by Tr/wrot, in effect letting the system advance through half a cycle of the field rotation, 
that also will not change the corresponding energy levels of the dressed state, which are after all defined over an entire 
period of the rotation. This transformation effectively connects measurements made for a > 0, above the mid-plane, 
to those with a < 0, below the mid-plane. The combined transform acts as follows: 
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The last of these is equivalent to n • ai, i.e., our symmetry operation would change the sign of the EDM energy shift. 
However, in the absence of this shift we can expect the following exact relations between the dressed state energies: 

(51) 

E{£,.ot, Brot, a, mp, ^) - E{£^ot, -B^ot, ~a, -mp, -fi) = 
E{£,ot,-B,ot,o^,mp,n) - E{£,.ot, B,ot, -a, -mp, -n) = 
E{£,ot,Brot,-(^,mp,n) - E{£rot, ~B,ot, a, -mp, - fl) = 
^{£rot,-Brot,-(^,mp,n) - E{£,ot, Brot, Oi, -mp, -n) = 0. 

Summing four equations and rearranging terms, we get that 

W''/\£,.ouB,oua) + W^/\£,ou-B,oua) + W^/\£,ou-B,ot,-a) + W^/\£,ouBrot,-a) = 0. (52) 

If we assume that the axial confinement is symmetric (not necessarily harmonic), and that our spectroscopy averages 
over an ensemble of ions oscillating in the axial motion with no preferred initial phase of the axial motion (we will 
later explore the consequences of relaxing this assumption) then the ions will spend the same amount of time on 
average at any given positive value of a as they do at the corresponding negative value of a, and thus the averaged 
results yield: 

< W^/\£,,u B.ot) > + < W^/'{£,ot,-B,,t) >= (53) 

The combined result, in the absence of exotic particle physics, is zero by symmetry. We did not need to invoke the 
various approximations that went into Eq. 1471 In particular, this null result is, unlike the traditional Berry's phase 
result, not based on the assumption of very small (i^rot/^mf^^rot)- Also, for conceptual simplicity we have discussed 
the result as being based on an average over quasi-static values of a, but the symmetry argument does not hinge on 
the axial frequency being infinitely slow compared to Wrot- 



E. Frequency- or Phase-Modulation of Axial Oscillation 



The trapped ions will oscillate in the axial direction at a frequency cj^, confined by an approximately harmonic 
axial trapping potential Uz = (l/2)Mcj|z^. Upon moving away from the mid-plane z — 0, the ions will experi- 
ence an oscillating axial electric field £z{t) = —Mu!'^z{t)/e. The geometric phase correction to the energy is then 
— mi^ojiotamax cos(a;2i), where amax = ^z.max/'S'rot is the maximum excursion of the tilt angle. Because the product 
WrotcUmax is again an energy, it is convenient to redefine the geometric energy contribution in terms of a frequency 

^max; 

Egco = WniaxC0s(w2i), (54) 

with Wmax = -SmpUJ^otama.^. 

For ujz — ^TT xl kHz, an ion cloud temperature of 15 K, an ion whose axial energy E^ is twice the thermal value, for 
Wiot and £rot as shown in the Appendix, then a transition such as W^^^ with Sm — 3 will have a maximum frequency 
modulation Wmax — 2t: x 400 Hz. 

Thus the electric field at the ion's location undergoes two motions, the comparatively fast radial rotation, and the 
comparatively slow axial wobble. We exploit the different time scales to create, for each instantaneous value of a, the 
rotation-dressed states worked out in the previous section. The effect of the axial wobble is then described by the 
time variation of the amplitudes in these dressed states. The time-dependent Schrodinger equation of motion for this 
is 




^gpHsBrot + cos{uj^t) ^ 

# -hpf^sBrot - COs(Wzt) 




(55) 



where a and b are the probability amplitudes for being in the \a) and \b) states, respectively. For typical experimental 
values, is about 27r x 1 kHz, Wmax will range as high as 27r x 1 kHz, and A (given by Eq. is perhaps 27r x 2 Hz, 
and S^BgpBiot is about 27r x 8 Hz. 

Eq. [55] describes a system again governed by a periodic Hamiltonian, and we will therefore follow a similar course 
to Sec. lIVDl and search for dressed-state solutions VP such that ^'(i -I- nr) = e^™'^^'(t). Of course, this will only be 
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FIG. 9: The apparent energy shifts between mp ~ +3/2 and mp = —3/2 states in upper (a,b) and lower (c,d) A-doublet levels 
versus Srot , "dressed" first by the electric field rotation (cjrot) and then by the ion's axial trap oscillation (iOz). At Brot = 0, 
there is an avoided crossing that mixes mp = ±3/2 states, with an energy splitting at the crossing of A"g^. In the limit 
SmpgFi-iBBrot 2> Acff , the dressed states are of good mp with an energy splitting slightly modified by Aefi. 



valid in the limit that Wiot S> Wz, a necessary condition to write the time-dependent Hamiltonian in Eg. 1551 First, we 
get rid of fast time-dependence by guessing solutions: 

a(t) = A{t) £ J„ (^) e--- (56) 



n— — oo 



n— — oo 



where are Bessel's functions of the first kind and A(t) and B(t) are slowly varying functions. We then substitute 
our trial solutions into Eg. [55] and use the recurrence relation (2n/a;) Jn(a;) = Jn-i{x) + Jn+i{x). We multiply through 



n — — 00 n' — — oo 



by 2^ Jn' ( 1 e'"-" ' or 2^ Jn' ( ) e*"^" ' as appropriate. We then integrate over an axial time 



period, 2t: /uj^, and make the approximation that A(t) and B(t) are unchanged over this small time interval. This 
approximation should be good as long as oj^ 3> A and Wz 3> gpfJ-sl^rot- The integration then yields. 




ot 



2 



ot 




(57) 



with 

= f; (If) ^» (^) A - A {^) A. m 

n— — 00 \ / \ / \ / 

This results in dressed-state energies, now as a function of yBrot, and not a, as seen in Fig. HI This clearly shows 
the requirement of S^FMs^rot > Aog in order to keep \a) and \b) as the dressed states. This is true despite the fact 
that an ion will sample the avoided crossing in Fig. |5] during its axial oscillation in the trap, as Wmax ^ SgF^BSrot 
in our experiment. Aoff will have a maximum value of A at Wmax/wz = and will oscillate about zero according to 

>/o('^max/Wz). 

For finite LOmnx/^zi the dressed states from Eq. 1571 onlv appear stationary if measured at integer multiples of the 
axial trapping period, 2tt/uJz. Consider states |+) and |— ), symmetric and antisymmetric combinations of states \a) 
and |6), respectively. In the limit that Smpgp^BT^rot ^ ^effj an ion initially in state |-|-) will oscillate between |-|-) 
and |— ) at the precession frequency ojq = {{3gp^^ ij,bB^^^ + deSes)^ + ^eff)^ > when measured at integer multiples of the 
axial trapping frequency. However, if our EDM measurement is made after a non-integer number of axial oscillations, 
or if the ions have different axial frequencies in the trap, the |-|-) to |— ) oscillation will be frequency modulated at 
■^max- For the example parameters, the frequency-modulation index Wmax/wz is less than 1, and thus the spectral 
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power of transition is overwhelmingly at wq, the quantity which symmetry arguments above show is unaffected by 
Berry's phase. In an ensemble of ions which have a random distribution of initial axial motions, the sidebands on the 
transition average to zero, and won't pull the frequency of the measured central transition. If instead the process of 
loading ions into the trap has left the ions with an initial nonzero axial velocity or axial offset from trap center, the 
measured frequency can be systematically pulled from loq. 

We note that increasing £rot or decreasing Wrot reduces the value of Wmax and thus the frequency modulation index. 
On the other hand, these changes also would have the effect of increasing the energy Erot of the micromotion of 
the ions in the rotating fields. For harmonic axial confinement, we find the frequency modulation for a dm — 3 
transition obeys the following relation cumax/'^z — 3(E2/Eiot)^^^. Thus to keep the modulation safely under unity 
for a comfortable majority of an ensemble of ions with an average Ej, given by T^, one needs to choose operating 
parameters such that E^ot > SOfc^Tz- This inequality in turn places stringent requirements on the spatial uniformity 
of Srot- On a time-scale slow compared to l/wrot, ^rot acts like a sort of ponderomotive potential analogous to the 
effective confining potential in a Paul trap. If E^ot = SO/csT^, then a spatial inhomogeneity in frot of only 1.5% 
already gives rise to structure in the ponderomotive potential comparable to Tz . 

To summarize the effect of axial motion: in the limit SgpfJ-sBrot > A, ions prepared, for instance by optical pumping, 
in state \a) (or \b)) will remain in \a) (or \b)). The energy difference between dressed states which are predominantly 
either \a) or \b) will be slightly modified by the avoided crossing. But the important combined measurement described 
by Eq. (|47p will continue to yield zero for — 0, and the sensitivity of that combined measurement to a nonzero 
EDM will not be much affected as long as Wmax/w^ < 1. 



F. Structure of the Measurements. What Quantities Matter 

In the remainder of this section, we will look at the possible effects of various experimental imperfections on our 
measurement. 

The symmetry argument in Sec. II V 5] presupposes the ability to impose a perfect ",B-chop", i.e., to collect data 
with alternating measurements changing quite precisely only the sign of Srot- If not only the sign but the magnitude 
of the rotating magnetic field alternates, the situation is more complicated. There will likely be contributions to the 
rotating magnetic field that are not perfectly reversed in our S-chop, including displacement currents associated with 
sinusoidally charging the electrodes that create the rotating electric field. These effects can be quantified with a value 
'^rot'^'^j and to lowest order they would appear as a frequency offset in the chopped measurement: 

W^/'{£, S^ot + Bfjn + W-I\e, -S^ot + ^^7") = ^gTl^BBf::^ ± 24fcff . (59) 

This offset is very nearly the same for the upper and lower states, to the extent that gp ^ gp, i.e., to the extent 
that 6gF = l/2{gp — g^) ^ gp. The effect of the stray field is reduced by combining measurements from the upper 
and lower states, in the form of a "four- way chop" : 

W^{£,ou Brot + ^^oD + W^iCu-Brot + Bf:n 

W'i£,ouB,.,t + Bfjn + W'i£,.,,~B,.ot + B^l^n 

= USgpfisBf^^'"'' + MeScff- (60) 

It may prove to be advantageous to shim the ,B-chop by deliberately adding a non-chopped rotating magnetic field, 
S*™, and adjusting its value until experimentally we measure 

w^iBrot + <r + b:T) + w^{~B,ot + + b:T) = o. (6i) 

Then, a measurement in the lower A-doublet state gives 

w\Brot + B:ir + Bi^n + w\~b,,, + bT' + Bt'^n = -2(1 + 4)de£.s, m 

9f 

yielding a still more accurate value for AdeScS- 

What we care about most then are: (1) Things that perturb VF" and differently, in particular the quantity 6gp, 
but also the quantity Sr], to be defined and estimated in Sec. IIVHI and (2) to a lesser extent, we care about effects 
which affect W"(i3rot) + W^"(— ^rot) the same way as they affect W'^{Biot) + W^{~~Brot), because, to the extent that 
they lead to a measurement 

I^"(6rot) + Wi-Brot) = +2defoff + SsysU (63) 
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we can mistake a nonzero value for (5syst as an indicator for a nonzero value of S^ot'^''. In that case, the shimming 
procedure discussed above to remove B^^^^^ would lead to a combined result from the four-way chop of AdeSeS + 
{A5gF I gF)5syst- This is down by a relative factor of {SgF/gp) compared to the effects that differentially perturb 
versus W^, but they could still be troublesome. And (3) to a still lesser extent, we care about imperfections 
that perturb individual measurements such as W"{Brot), even if they do not perturb the S-chop measurement, 
W'^{Brot) — W^{—Brot), because, to the extent that they vary over time, or depend on the trajectory of an individual 
ion in the trap, they can reduce coherence times. This leads not to systematic errors, but to a reduction in the overall 
precision. 

In addition to the B chop, state chop, and four-way chop discussed above, we can perform a rotation chop, by 
changing the sign of Wrot • Our hope is to keep experimental imperfections to a level where the four- way chop is by 
itself already good enough to suppress systematic error below the desired level. Then repeating the entire series of 
measurements with the opposite sign of Wrot (rotating the field CW instead of CCW) will to the extent it yields the 
same final value of AScsde provide a useful redundant check. 



G. An Estimate of Sqf^ 



3/2 



There are two leading contributions to 5(?f=3/2 in our molecule. In the regime in which we will operate (a regime 
wherein f2 is a signed quantity) they are to a good approximation independent of each other. These two contributions 
are the zero-field difference and the induced difference caused by the applied electric field. In the zero-field limit, the 
former is dominant. However, in the limit in which we are working, the latter dominates. 

The zero electric field contribution arises due to centrifugal distortion effects in the molecular Hamiltonian. In 
Sec. IIVBI we wrote the Zeeman Hamiltonian in Eq. [211 We omitted two terms which connect states of il — il. 
Therefore, these states will give rise to a parity dependent g-factor for each J-level. The Hamiltonians which govern 
this interaction can be found with the us e of perturbation the ory i n a manner similar to the approach used to find 
the A-doubling parameters. Brown et al. [ll8j | and Nelis et al. [l20j | have written these terms as 



^^ZocmanDist — — -grS l^'B {B+J- S+ S - + B- J+ S - S+) 
-H^ZocmanDoub — -glgfJ.B{B+J+S^+B-J-S^), (64) 

where -ffzcomanDist IS the centrifugal distortion induced by the magnetic field and i?zccmanDoub is the Zeeman induced 
A-doubling. i?zoomanDist IS parity independent while -ffzccmanDoub is parity dependent. Due to the nature of the 
perturbation approach, we can estimate the size of g'^g in terms of the A-doubling J = 1 energy splitting Wof 



2B, 



(65) 



In addition, if the '^Ai state of interest is composed of a {s)a{d)S molecular orbital (where (s) and (d) refer to atomic 
orbitals with I = 0,2), then grs = g'rs is expected. The difference in zero- field g- factors is then given by twice the 
value in Eq. (|55|) . It is evident that this effect is quite small, of the order 10^^ for HfF+. 

The electric field dependent g-factor arises due to the mixing of rotational levels J in the molecule. The levels with 
J = 2, while far away in energy compared to the Stark energy dmffrot, are perturbers. In the signed basis, the mF 
sub-levels in the J = 2 level have a smaller •jf value than do the mp sub-levels in the J = 1 level. Therefore, the 
states which go up (down) in energy in the J = 1 level "gain" ( "run" ) on (from) the J = 2 level. When one includes 
the effects of Hyperfine interactions, there are multiple connections to each sub- level. In the J = 1, niF = ±3/2 levels 
that we are interested in, we can write an analytic expression for the electric field dependent SgF factor 

. . <fgrot ffF rr. r./ t Tn2 f F 1 F' V ( J I J' V j F' J' I V 

S9As.o.)^J:^^-^^^[f,f,j,j] oruF) [-no n) [j F ij ^ ^''^ 



where [J, J', . . .] — \J (2J 4- f )(2J' -I- 1) . . .. The terms in parentheses are 3 J-symbols while the term in curly brackets 
is a 6J-symbol. The sum runs on all states connected to |J, F) by the electric field. In the case of HfF+ with a 
J — \,F = 3/2 ground state, the sum contains the J' = 2 and F' = 2,12,^/2 states. Since the rotation constant 
Bf. is far larger than either dmfS'rot or Ehf, only is included in the perturbative expression for 5gF{£mt)- For the 
parameters here, this contribution is 

^5^=3/2,^ X Qdmffrot 

^('^rot) = , (67) 
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which means that the fractional shift 5gp^^/2/ gF=z/2 is a few 10~^. The same approach gives that the electric field 
"g" factor, 7^, will shift in the same manner such that 5^f/if ~ 10^'*. 

For rotating fields, another contribution to 5gF arises from non-vanishing value of i^rot/ [dmi^rot)- The states 
with J7 = 1 and f2 = — 1 are equally affected by the rotating field since they have an equal Stark shift in the 
absence of hyperfine interactions. However, because the levels with \F = 3/2, mi;- = ±1/2) are repelled by the lower 
\F — 1/2, mi? = ±1/2) states, the effective Stark difference between to_f levels with 51 = —1 (upper levels) is smaller 
than the same niF levels with VI — +1 (lower levels). The scale at which this difference will appear is then determined 
by how much the lower hyperfine state pushes on the upper due to the coupling induced by the electric field. 

'-^^ f , f\ . (68) 

This fractional shift is of the order a few 10~^ and is therefore about the same magnitude as the electric field induced 
mixing of higher rotational levels. 



H. Dependencies on Srot 

Proximity to the avoided crossing shown in Fig. [5Jb) means that the transitions and will have residual 
dependencies on frot, which in turn may lead to decoherence or systematic errors. We characterize the sensitivity of 
W"/^ to small changes in E-not with the following expansion 

W^"/^(Ct + (Jfrot, Srot) = T4^"/^(£:°ot, ^rot) + TI^'^SE^ot (69) 

with 



5f MB^rot'S'r 



(70) 
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using the expressions in Eqs. (j42|) and (|46p. Any spatial inhomogeneity in frot that does not average away with ion 
motion will lead to a decoherence rate given approximately by r/S£mt- 

In terms of systematic errors, if chopping the sign of Brot gives rise to an unintended systematic change in the 
magnitude of S^ot (call it (5fchop), for instance due to motional fields discussed later, or due to ohmic voltages generated 
by the eddy currents, then there will be a frequency shift in a ,B-chop combination, 2ri'^/^S£chop- To the extent that 
Sr] = ^(?/" — Tj^) is nonzero, some of this shift will survive a four- way chop as well. The dominant contribution to St] 
is likely from 6A, rather than from 6gF- Assuming this limit, the systematic error surviving is 

8^rjS£chop « 6 ^-^^^^^ Jy'J'S^chop- (71) 

For a large but not inconceivable value for (5£chop of 100 /iV/cm, and for other values as in the Appendix, this works 
out to comfortably less than 100 /xHz, and is therefore not a problem. But this error would scale as £^^, and thus 
could cause trouble if for other reasons we chose to decrease frot ■ The science signal is roughly independent of E^ot , 
which should allow for the source of error to be readily identified. 



I. Perpendicular ^?-Fields 

The quantization axis is essentially defined by fiot- The shift of the various levels |a), \b), |c), \d) due to a component 
of the magnetic field perpendicular to £i-ot is on the order of 

^ 3 jgF^iBB^)' .^2) 

4 jFdmiSrot 

for the upper/lower states. In the absence of rotation, the lowest-order correction to iy'/^(Srot) goes as 

_ 3 glii%BlB,ot 

2 (7FC?mffrot)^ ' 



(73) 
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For reasonable experimental parameters, this will be a negligible number. The lowest-order correction to the state-chop 
combination, VF"(S,.ot) - W^iB, ot) is smaller still and goes as 

3ff|4^i^^ ^^^^ 

It is similar in form to the difference in (/-factors caused by the rotation of the field. 

When we turn on rotation, there is an additional larger contribution to W/^(Srot)- If we assume (as a worst case) 
that B± is purely axial, not azimuthal, then the lowest-order effect of B± is to tilt the quantization axis by angle 
given by 

±gFi^BB± ^^^^ 



with the +(-) in the numerator corresponding to the upper(lower) states and the -l-(-) in the denominator corresponding 
to the = -I(+l) states. This has the leading order effect on 14^"/^ of 

^i^rotgFfJ.%B±Brot ^^g^ 
(7Frfmffrot)^ 

even a rudimentary nulling of the Earth's magnetic field, say to below 25 mG, will leave this term negligible, for 
parameters in the Appendix. Its contribution to the state chop, VK"(Srot) — W^{Bi-ot), is still smaller by dmti^^rot /Ehf 

(^rotglfJ.%B±Brot ^^^^ 

7FC?mffrotEhf 



J. Stray Contributions to : Uniform or Time- Varying B Fields 

In the previous section we have seen that the effects of Bj_ are small. Spatial or shot-to-shot variation in i3| | , on the 
other hand, can limit coherence time through its contribution to W^l^ . The biggest contribution to B\\ is of course 
the intentionally applied rotating field B-^ot. Let's examine the various other contributions to B\\. 

Static, uniform fields: B fields of this nature are relatively harmless. B\\ is defined relative to the quantization 

axis £iot- The time-average of B\\^ is (S • £rot)- Since Sj-ot sweeps out a circle with angular velocity Wjot, the 

contribution to the time-averaged B\\ from a uniform, static magnetic field averages nearly to zero in a single rotation 
of the bias electric field, and still more accurately after a few cycles of axial and radial motion in the trap. The average 
electric field in the ion trap must be very close to zero, or the ions would not remain trapped. In the case of certain 
anharmonicities in the trapping potential, however, one can find that the average value of f rot is nonzero, even if the 
average value of frot is zero. For instance, an electrostatic potential term proportional to z'^, along with a uniform 
axial magnetic field B^ , will for an ion with nonzero axial secular motion, yield a nonzero (S| | ) • In addition, nonzero 

Bz will interact with the tilt of £rot oscillating with an ion's axial motion at to cause a frequency modulation 
similar to the one discussed in Sec. lIVEl A uniform magnetic field in the x-y plane will cause a frequency modulation 
at Wrot- If the modulation index for either of these modulations approaches one, the modulation will begin to suppress 
the contrast of spectroscopy performed at the carrier frequency. For uniform magnetic fields with amplitude less 
than 10 mG (achievable for instance by roughly nulling the earth's field with Helmholtz coils), frequency modulation 
indices will be small, and, barring pathologically large electrostatic terms, the mean shifts from uniform, static B 
fields will be less than 1 Hz and can be can be dealt with by means of an applied /B^o™ as discussed in Sec. IIVFI 

Time- varying magnetic fields with frequency near Wrot can cause more trouble. If the time between the two Ramsey 
pulses used to interrogate the frequency is iRamsoy, then the dangerous bandwidth is 1/iRamsey, centered on Wrot- We 
discuss in order (i) thermally generated fields from the electrodes, (ii) ambient magnetic field noise in laboratory, (iii) 
magnetic fields associated with the application of frot, oscillating coherently with £rot, (iv) shot to-shot variation in 
magnitude of applied Srot, and (v) spatial inhomogeneities in Srot- 

(i) Proposed EDM experiments on trapped atomic species such as Cesium are vulnerable to m agnetic field noise 
generated by thermally excited currents in conductors located close to the trapped species [l26l |. In our case, the 
effect is less worrisome because, vis-a-vis the trapped atom experiments, our bandwidth of vulnerability is centered 
at much higher frequency fields, because our molecules are trapped considerably further from the nearest conductors, 
and because the sensitivity of our measurement of de to magnetic field noise, which goes as gFlJ^Bl^'e'R is d own by 
a factor of 10*. The spectral density of thermal magnetic field noise (which is calculated in reference jl27l | in the 
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simplified geometry of a semi- infinite planar conductor) will surely be less than 1 pG/Hz^/^ in our bandwidth of 
vulnerability. This effect is negligible. 

(ii) Like thermal magnetic noise, technical magnetic noise in our lab arising for instance from various nearby 
equipment will not so much decohere an individual measurement as generate shot-to-shot irreproducibility between 
measurements. What level of noise are we sensitive to? As we discuss in Sec. IVII below, the precision of a single trap 
load is unlikely to be better than 300 mHz, meaning magnetic field noise less than 0.2 //G/Hz^/^ won't hurt us, for a 
1 s interrogation time. Measurements made in our lab show that there are a number of magnetic field "tones" of very 
narrow bandwidth, associated with harmonics of 60 Hz power and various power supplies. As long as we choose Wrot 
to not coincide with one of these frequencies, in the range of 50 kHz to 300 kHz ambient magnetic frequency noise 
in our lab has spectral density typically less than 0.02 fiG/Hz^^^. For this reason, at least for the first generation 
experiment, there will be no explicit effort to shield ambient magnetic field other than to use Helmholtz coils to 
roughly null the earth's dc field. The steel vacuum chamber will in addition provide some shielding at 100 kHz. 

(iii) In traditional eEDM experiments, one of the most difficult unwanted effects to characterize and bring under 
control is magnetic fields generated by leakage currents associated with the high voltages on the electrodes that 
generate the principal electric field. In our case the bulk of the electric field £cS is generated inside the molecule. The 
laboratory electric fields are measured in V/cm, not kV/cm, and leakage currents as traditionally conceived will not 
be a problem for us. On the other hand, the electric field does rotate rapidly, and thus the electrode potentials must 
constantly oscillate. Displacement currents in the trapping volume between the electrodes, and real currents in the 
electrodes themselves and in the wire leads leading to them, will generate magnetic fields with spatial gradients and 
strengths that oscillate coherently with frot at the frequency ojiot- 

The spatial structure of the oscillating magnetic fields will depend on the geometry of the electrodes and in particular 
on the layout of the wire leads that provide the current to charge them. In principle, shim coils can be constructed just 
outside the trap electrodes and driven with various phases and amplitude of current oscillating at oji-ot, all in order to 
further control the shape of the magnetic field. The one immutable fact is the Maxwell equation, V x S = c~^d£ /dt. 

The dominant time dependence of the electric field is from the spatially uniform rotating field, and thus for a 
circular field trajectory, the dominant contribution to the magnetic field structure goes as 

VxB = ky' (78) 

with 

k ^ c-H„,.„, ^ 350 „G/c,„ X (^) , (79) 

where y' is the direction in the x-y plane orthogonal to the instantaneous electric field. 

The curl determines only the spatial derivatives of B: B itself only depends on the boundary conditions. An idealized 
arrangement of current carrying leads and shim coils could in principle force the B field to be 

6'^'=^' = kx'z. (80) 

where k is given by Eq. [79] and x' is displacement in the x-y plane along the direction of the instantaneous rotating 
electric field. These fields would be perpendicular to the quantization axis provided by the electric field, and would 
have negligible effect on the transitions of interest. 

While realizing such an idealized displacement field would be very difficult, there are relatively simple steps to 
take to minimize the displacement fields. For instance, each rod-like electrode can be charged up by two leads, one 
connected to each end of the rod, with the leads running along respective paths symmetric in reflection in the z=0 
plane to a common oscillating voltage source outside of the vacuum can, at z=0. It is worth considering a maximally 
bad electrode layout, to put a limit on worst-case performance. Our electrodes will be spaced by about 10 cm and 
mounted in such a way that their capacitance to each other or to ground will be at worst 5 pF. If the charging current 
is provided entirely by a single lead connected to one end of the rod, the peak current running down the rod near its 
center will be 80 /xA, leading to a worst-case field magnitude at the trap center of about 20 /iG, and a contribution 
to W^/^ of perhaps 2.5 Hz. Spatial gradients of this effect, and shot-to-shot irreproducibility of this effect will not 
contribute to decoherence at the 0.1 Hz level. As for its contribution to systematic error, this shift will survive the 
S-chop, but will be suppressed in the four-way chop by the factor (Sgp/gF), perhaps a factor of a thousand. For still 
better accuracy the shift should be nulled out of the B chop by adjusting ySf^™, as discussed in Sec. lIVFl 

(iv) Given that the main effect of B^ot is to apply an offset frequency, igpl^BBmt of perhaps 8 Hz, and given that 
(see Sec. IVip the single-shot precision is unlikely to be any better than 300 mHz, the shot-to-shot reproducibility of 
Slot need be no better than a part in 30, a very modest requirement on stability. Decoherence then is not a problem, 
but a potential source of systematic error arises if the the B chop is not "clean" that is if Srot before the chop is not 
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exactly equal to — Srot after the chop. This sort of error could arise for instance from certain offset errors in op-amps 
generating the oscillating current. Experimentally, one adjusts B^^l™ to cancel these offsets, but even in the absence 
of that procedure, the four-way chop cleans up these sorts of errors. For a rather egregious fractional deviation from 
;B-chop cleanliness of, for instance, 1%, and for {Sgp/gp) < 0.001, the systematic error remaining after the four-way 
chop is 10^^ of the offset frequency of perhaps 8 Hz. In HfF+ this is a systematic error on de of 10~^^ e cm. For 
ThF+ the error as referred to de is smaller still, and of course if we avail ourselves of ^B^Jj™ so as to null the post 
S-chop signal to <100 mHz, the systematic error on de will be less than 10~^^ e cm for either species. 

(v) The largest single contribution to decoherence (with the exception of spontaneous decay of the '^Ai line to a 
lower electronic state) will likely be due to spatial inhomogeneity in the applied rotating bias field B^ot ■ That is to say, 
spatial inhomogeneities in B that rotate in the x-y plane at frequency Wrot- First-order spatial gradients in Biot are 
not important, because ion secular motion in the trap will average away the effects of these gradients leaving only the 
value of Biot at the center of the trap. Second-order spatial gradients on the other hand will lead to nonzero average 
frequency shifts whose value will vary from ion to ion in a thermal sample of ions, depending on conserved quantities of 
individual ion motion like the axial secular energy or radial secular energy Ep, quantities with thermally averaged 
values of kT^ and kTp, respectively, and with ion-to-ion variation comparable to their mean values. The B^t will be 
generated by current-carrying rods which are of necessity within the vacuum chamber because of the screening effects 
of a metal vacuum chamber. Unless particular care is taken in the design of these rods, the second-order spatial 
gradients in Brot will scale as where X is the characteristic size (and spacing) of the current carrying rods. 

The contribution to the inhomogeneity of the time-averaged value of B^ot experienced by a thermal sample of ions 
orbiting in a cloud with r.m.s size r is then of order (r'^ / X'^)Brot, leading to an ion-to-ion frequency variability of order 
{r^ / X'^)ig F HBB^ot For planned parameters of the experiment, (r^/X^) is of order 0.01. We have seen from Sec. lIVEl 
above that the quantity SgpfJ-sBrnt must be at least about five times larger than A in order to make the eigenstates 
in the rotating fields be states of good niF. Thus in the absence of explicit apparatus design to null the second-order 
spatial gradient in Srot (The rod-like electrodes that bear the charge that generates frot are in the second-generation 
trap the same objects that carry the current that generates Srot and thus their shape is already subject to multiple 
design constraints) we may have to live with a decoherence rate from this effect on the order of 0.05A, perhaps 0.5 
s~^, for the experimental values given in the Appendix. 

The inhomogeneity in B^ot should reverse quite cleanly with the B chop, and residual imperfections there will be 
cleaned up with the four-way chop, and thus the effects of the second-order gradients in Srot are expected to be 
predominantly a source of decoherence, rather than systematic error on measured dg. 



K. Stray Contributions to Static B-Field Gradients 

We now return to discussing static magnetic fields, now including the effects of spatial gradients. With the charac- 
teristic size of the ion cloud r being smaller than the characteristic distance X from cloud center to source of magnetic 
field by a ratio of 0.1 or smaller, it makes sense to expand the field about the uniform value at the trap center. The 
most general first-order correction to a static magnetic field in the absence of local sources can be characterized by 
five linearly independent components as follows: 

^ = ^axgrad(^^ ~ f ^ " 2^^ 

+Sira„s(2;i - yy) 

+B'i(yx + xy) 
+B'2{zx + xz) 

+B',iyz + zy) (81) 

By far the most important effect of these terms is the "micromotion-axial gradient interaction." As discussed in 
Sec. IIII B] above, the displacement of an ion's circular micromotion r^ot is exactly out of phase with the rotation of its 
quantization axis £, see Eq. [HI Averaged over a cycle of Wrot, this will give rise to a nonzero average contribution to 
B\\ and cause a shift in W"^^ given by 3gFM-B'S^xgrad'"iot — 3gFMS'BaxgradS^i-ot/(Afw^ot). A guess for a possible value 
of stray B'^^^^^^ is 2 mG/cm, which for anticipated experimental parameters would lead to a shift in W^^^ of order 4 
Hz, and this shift would survive the B chop. As with the effect of displacement currents, one expects the systematic 
effect of the shift to be reduced after the four- way chop by (Sgp/gp), but for maximum accuracy the effect should 
be shimmed out of the B chop, either by adjusting the value of B^^l™, or by applying (say with anti-Helmholtz coils 
external to the vacuum chamber) a compensating value of ^axgrad- 

A smaller effect arises from the interaction of the magnetic field gradient with the component of the electric fields 
responsible for providing ion confinement, which after averaging over cycles of Wrot and Wjf, always point inward, 
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giving rise to a net inward-pointing time average oi £ . If we look at only the component of the first-order magnetic 
field gradient that points towards or away from the trap center 

^central ^ (g,^^^^^ _ B'.,^^,,J 2)xx + {-B[,^^, - B',^^,.,j2)yy + B',^^,^^zz. (82) 

The net contribution to B\\ comes from integrating, along the rf and rotation micromotion trajectories, over first a 
rotational cycle, and then an rf cycle, and then a secular cycle in a given direction. We assume that the trap is suffi- 
ciently harmonic that there is no cross-dimensional mixing of secular energy, that ljx, Wy, and uiz are incommensurate 
and with principle axes as defined in Eq. 1811 and that hard-momentum-changing collisions are rare enough so that, 
during the duration of a spectroscopic measurement, there is no change in E^, the sum of the kinetic and potential 
energy associated with an individual ion's secular motion in the «th direction. The contribution to B\\ is then, 

5B\\ = -{{B[,,J2 - B;,g,,d/4)E. -t- {-B[,^J2 - BLgrad/4)E, -I- 6Lg,,dE.)/(ef,.ot). (83) 

The contribution to B\\ averaged over a thermal sample of ions is given by the above expression with E,; replaced 
by Ti. Note that for = Ty — Tz, several terms cancel and the thermally averaged contribution to B\\ is just 
'^axgradEz/lSe^rot)- The dccohcring effect is comparable because within a thermal sample, E^:, Ej,, and E^ will in 
general differ from one another for a given ion, and between different ions. For ;Baxgrad ^"^^ '^trans sach about 2 
mG/cm, ion temperatures about 15 K, the mean shift in W"/^ for typical experimental parameters given in the 
Appendix might be 30 mHz, with a comparable contribution to dephasing. 

The three remaining terms in the first-order gradient, B'l, B2 and B';^, will contribute to a shift in B^ only when 
combined with other (usually small) trap imperfections, for instance the plane of rotation of S^ot being tilted with 
respect to the principal axes of the confining potential. The net effects will be correspondingly smaller than those 
from 

Just as with the second spatial derivative of B^ot , the spatial derivative of ^axgrad ' coupled to a thermal spread in 
the size of ion orbits, can give rise to decoherence. Of course, Kaxgrad is defined already as a first spatial derivative 
of a magnetic field, thus the dephasing arises from a third derivative of the field, and its rate should be down from 
the mean size of the shift (roughly estimated above at 4 Hz) by a factor of order {r/X)'^, or a factor of one hundred. 
Even spatially uniform -B^xgrad could give rise to decoherence if there is a spatial dependence in Trot- The fractional 
change in rrot is the same as the fractional change in fj-ot- As discussed in Sec. lIIIBl this should be smaller than 0.5% 
over the typical size of the ion sample. 

As a coda to this subsection, it is worth considering that applying a very spatially uniform Biot may be very 
challenging because of difficult-to-model eddy currents induced in electrodes and light-gathering mirrors. On the other 
hand a purposely applied static ^a^giad would be perturbed only by the magnetic permeability of trap construction 
materials, which can be minimized and modeled. One way or another we will need to bias away from the avoided 
crossing discussed in Sec. lIVDi but it may turn out that this can be accomplished with greater spatial uniformity 
and thus with a lower total decoherence rate by omitting the applied B^ot altogether, and providing the bias with a 
deliberately applied ^axgrad field. The B chop could be accomplished by chopping the sign of B'^^^^^^. The parity 
invariance argument of Sec. IIVDI above can readily be modified to describe a chop of ^Saxgrad rather than a chop in 
Biot- 

To sum up subsections IV. J and IV. K, we have looked at a range of ways in which various contributions to B\\ can 
shift W" and W^. Decoherence due to shot-to-shot fiuctuations or spatial inhomogeneity should not be a problem 
out to beyond 1 s coherence times. Various effects can shift W" and by as much as a few Hz, and this shift can 
survive a B chop. With Sgp/gp on order of 10"'', and See estimated at 90 GV/cm in ThF+, after a four-way chop 
the remaining systematic error will be a few 10~^^ e cm, but this can be dramatically reduced by tuning away the 
post-S-chop signal with B^^^'^. The most dangerous systematic error would be if Byot were systematically different 
between measurements on the upper and on the lower states. Chopping between upper and lower states will be 
determined by variations in optical pumping, which should be well decoupled from the mechanisms that generate 

Brot- 



L. Relativistic (lon-Motion-Induced) Fields 

The largest component of the velocity on the ions is that of the micromotion induced by Siot] for reasonable 
experimental parameters it will be less than 1000 m/s. In typical lab-frame magnetic fields of a few mG, the motion 
will give rise, through relativistic transformation, to electric fields of order of a few ^V/cm, which are irrelevant 
to our measurement. Conversely, motion at 1000 m/s in typical lab-frame electric fields of 10 V/cm generates a 
magnetic field of 0.1 /iG. This field will be rigorously perpendicular to the electric field, the quantization axis, and 
thus represents only a negligible modification to the generally unimportant B±. 
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FIG. 10: Over one rotation of frot, both S^i and £z are quasistatic. The total electric field is the sum of all three and its 
trajectory over one cycle of £rot is plotted as the dotted line projected onto (a) the x-y and (b) the x-z planes. The electric 
field trajectory is a circle of radius £rot, parallel to and displaced from x-y plane, a circle whose center is offset from the z-axis 
by frt. In the limit | £^rf | ^ | £rot |, the solid angle subtended from the origin by this circle differs only slighlty from that 
subtended by a circle with vanishing £,■{. The magnitudes of both £,■{ and £z relative to £iot are very much exaggerated for 
clarity. 



M. Effect of RF Fields 



The effects of the rf electric fields providing Paul trap confinement are best understood by putting them in the 
context of a three-tier hierarchy of electric field magnitudes and frequencies. 

(i) frot, the nominally uniform, rotating electric field, with field magnitude of perhaps 5 V/cm and frequency 
cjj-ot = 27r X 100 kHz. 

(ii) £rf, the Paul-trap fields, are highly inhomogeneous, but at a typical displacement in the x-y plane of perhaps 0.5 
cm, the field strength might be 75 mV/cm, or two orders of magnitude less than that of S^ot, oscillating at a frequency, 
Wrf = 27r X 15 kHz which is one order of magnitude less than uj^ot- At a fixed point in space, the rf fields average 
rigorously to zero over time, but averaged instead along an ion's rf micromotion trajectory, the rf fields contribute to 

(iii) the inward-pointing trapping electric field, again very inhomogeneous but with typical strength down from peak 
rf-field values by factor of (wrf/w^), another order of magnitude, to perhaps 5 mV/cm. From the ion's perspective, 
the direction of the trapping fields oscillate with the ion's secular motions in the trap, at frequencies iOi of perhaps 
27r X 1 kHz, the slowest time scale by an order of magnitude. 

The effects of the strong, fast f rot have been discussed extensively throughout Sec. IIVI and those of the weak, slow 
trapping fields were covered in Sec. IIV Kl above. In this subsection we argue that the rf electric fields, intermediate in 
both frequency and strength, are the least significant of the three categories. 

The effects of the rf fields averaged over the rf micromotion trajectory are discussed in Sec. IIVKI The remaining 
part averages to zero in one rf cycle, but is roughly frozen at a single value over the duration of one cycle of uj^ot- The 
dominant source of the rf fields' time-averaged contribution to transitions W"/' is in very small corrections to Berry's 
phase energy associated with the rotation of f rot- See Fig.[TUl The correction to the solid angle arising from Erf goes as 
(E2/Erot)(Erf/Erot)^- If wc iucludc a factor of Wrot to get a Berry's energy shift and evaluate for typical experimental 
parameters, the magnitude of the resulting frequency shift will be about 20 mHz, and will oscillate in sign with the 
axial secular motion. The magnitude of radial rf fields scales linearly with the radial secular displacement. If secular 
freqencies were commensurate, in particular if ojz = 2ujr, then this 20 mHz shift could contribute to a decoherence 
rate at the negligible level of a few tens of mHz. For incommensurate ratios of Wz/w^; or Wz/wy, the rf fields will be 
still less important. 



N. Systematic Errors Associated with Trap Asymmetries 

The symmetry argument of Sec. IIV Dl was based on parity invariance. This argument is only as good as refiection 
symmetry of the electric and magnetic fields in the region of the trapped ions. In this section we look, as an example, 
at the consequences of a symmetry imperfection. 

The electrodes used to generate £rot have been numerically designed to make frot as spatially uniform as possible, 
but imperfections in design and construction of the trap and imperfect drive electronics will lead to some residual field 
nonuniformity. Suppose that the magnitude of the f rot was consistently larger in the region of the trap for which z>0. 
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so that the value of £^ot over the z>0 half of an axial secular oscillation is about 0.3% larger than that experienced 
over the z<0 half. Thus the frequency modulation of perhaps ± 500 Hz, discussed in Sec. lIVEl will no longer average 
to precisely zero over an axial cycle but instead a net contribution of about 1.5 Hz to W". Such a frequency shift 
would survive a B chop, and, following the protocol discussed in section IV. F, we could very likely incorrectly identify 
this shift as arising from the presence of a B^*^'*'*', and apply a value of B**'"'" to largely null the 1.5 Hz shift. After 
a complete four-way chop, we would be left with a systematic error on the order of [Sgp / 9f)'x^-^ Hz, or about 0.4 
mHz. 

For the value of £cs estimated for HfF+, a 0.4 mHz error corresponds to a systematic error on de of the order 
of a few 10~^^ e cm. For ThF+, the error on de would be about three times smaller. We continue a more general 
discussion on systematic errors in Sec. IVI D] below. 



V. COLLISIONS 



The overarching strategy of the trapped-ion approach to precision spectroscopy is to accept low count rates in 
exchange for very long coherence times. In some previous precision measurement experiments with trapp ed ions, the 
very best results ha ve co me from taking this to the extreme limit of working with only one ion 128l4l35l | , or in some 



cases a pair of ions [136j, in the trap at any given time. More often however, optimal precision is achieved working 
with a small cloud of trapped ions. In this section we evaluate various detrimental effects of ion-ion interactions. 



A. Mean-Field 



With no electrons present to neutralize overall charge, even a relatively low density cloud of ions can have a 
significant mean-field potential. A spherically symmetric sample of iVion ions confined within a sphere of radius r will 
give rise to a mean-field potential 

f/mcan-ficld - , 3 ^ x f ^\ ( -^—) . (84) 



kB VlOOOy V0.5 cm 

At values of the mean-field interaction energy comparable to or larger than fc^T, there is a risk of instabilities, viscous 
heating, and other undesirable effects; even in their absence, systematic errors are more difficult to analyze in the 
strong mean-field limit. Ion-trap experiments have been performed at much higher mean-field strengths, and indeed 
there have been precision spectroscopy experiments done in systems for which the interaction potential even between 
an individual pair of nearest- neighbor ions is much larger than ksT . However, these systems exhibit a high degree of 
spontaneous symmetry breaking including crystallization jl37| . 

For the purpose of this paper, we assume the experiments will be done in the low mean-field limit, say 

C^mcan-ficld ^ \^bT. (85) 

In this limit, mean- field effects are relatively benign, and can be modeled as a modest decrease in the trap confining 
frequencies, Wi, plus the addition of some anharmonic terms to the potential. Crucially for the arguments presented 
in Sec. IIV D[ these additional modifications do not break any of the reflection- or rotation-based symmetries of the 
trapping fields. We note that Egs. IMl and 155] combine to set limits on various combinations of the ion number, A^ion, 

ion temperature, T, cloud radii, cx ^J^^^ and mean ion density, n oc . This necessitates making various 

compromises in selecting operating parameters. 

In Sec. IIV El we saw that the axial component of the electric field at the ion's location, f^, tilts the rotating electric 
field and gives rise to an apparent shift of the energy of our spectroscopic transition, linear in E^- This energy shift 
integrated over time in turn gives rise to an oscillatory phase shift, A0 = J S^dt. In a one-component ion cloud, 
the effects of long-range, grazing-angle ion-ion collisions may be thought of as simply a fluctuating component to the 
local electric fleld, and the integrated effect of those fluctuations will make a random contribution to the phase shift. 
We present a simple argument to show that the resulting rms spread in phase does not continue to increase with 
time but reaches a steady-state asymptote. This is because £z not only shifts the transition energy, it also causes 
an axial force and corresponding acceleration, which, like the shift in transition energy, is linear in S^. Integrated 
over time, Ap^ = e J £zdt, this fluctuating force results in a fluctuating momentum. But we know that the combined 
effect of a trapping fleld and a large number of random collisions will not cause the rms momentum to randomly walk 
without bound but rather to be loosely bounded by a characteristic thermal value, \/(pI) ~ -s/Mfc^T^. This is the 
nature of the thermal equilibration process - once an ion has developed a super-thermal momentum, further collisions 
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FIG. 11: Geometric phases accumulated during an ion-ion collision, (a) A typical ion-ion collision trajectory (red), resultant 
Rutherford scattering angle, 6, and ion-ion interaction electric field, fion, are shown in the coUision plane (blue). For clarity, the 
colhsion plane has been taken perpendicular to the instantaneous direction of frot. (b) During an ion-ion collision the molecular 
axis adiabatically follows the net electric field vector, £rot -l-fion, and traces out the contour (black) on the unit sphere (yellow). 
The solid angle, A^(6), subtended by this contour gives rise to a geometric phase accumulated by the eigenstates during the 
collision. This leads to decoherence of the spectroscopic transition, see text. 



are biased to reduce the momentum. Since both the phase excursion and the momentum excursion are linear in the 

time-integrated axial electric field, we can estimate -y/^A^^ ^^\/ (Pz) ~ \J 2E^l ' ^S^i'^' discussed in Sec. lIVEl 

if Eiot ^ SOksTz, the phase fluctuations for each ion's spectroscopic transition will be bounded by a value less than 
one radian, so that there will be no loss in spectroscopic contrast in a Ramsey- type experiment. 

The argument in the paragraph above hinges on the assumption that the electric field arising from the ion cloud's 
mean- field distribution and from grazing-angle collisions is small in magnitude compared to Srot, so that the shift in 
Berry energy is linear in the axial component of the electric field. For higher values of the ion temperature or lower 
values of £rot, a pair of colliding ions can get so close to each other that the electric field is, transiently, comparable 
to or larger than £^ot- We discuss the consequences in the next subsection. 

B. Geometric Phases Accumulated During an Ion-Ion Collision 

As discussed in Sec. lIVDl when a spin adiabatically follows a time- varying quantization axis it acquires a geometric 
(Berry's) phase. For the eigenstates in Fig. [Ub), the geometric phase factor can be written as exp {zLimpA), where 
A is the solid angle subtended by the contour on the unit sphere traced out by the time-varying quantization axis. 
Thus, the relative phase generated between the \F — 3/2, mp — ±3/2) states used for spectroscopy is (j) = 3A. The 
concern of this subsection is how ion-ion collisions cause uncontrolled excursions of the quantization axis leading to 
random geometric phase shifts and decoherence between spin states. These uncontrolled phase shifts will be written 
as A(f) — 3AA to distinguish them from the calibrated geometric phases in the experiment. 

The instantaneous quantization axis for the molecular ion eigensates is defined by the net electric field vector at 
the location of the ion. During a collision, this axis is defined by the vector sum of the rotating electric field, £i-ot, and 
the ion-ion interaction electric field, fion- Both of these are time-varying vectors, however typical ion-ion collisions 
have a duration short compared to the rotation period of f rot so for the purpose of this discussion frot will be taken 
as stationary. Thus, the problem is reduced to calculating the excursion of the quantization axis under the time 
variation of Sion- A typical ion-ion collision is shown in Fig. Illf a) and the effect of this collision on the quantization 
axis is shown in Fig. [TITb). 

At the temperatures of our trapped ion samples, no two ions are ever close enough for the details of the intermolecular 
potential to matter. Only monopole-monopole and monopole-dipole interactions matter. Further, the translational 
degree of freedom may be treated as purely classical motion in a 1/r ion-ion potential, with the initial condition of 
a given coUisional event characterized by an impact parameter and relative velocity. The outcome of the collision 
depends not only on the magnitudes of the impact parameter and of the velocity, but also on their angles with respect 
to the ambient electric bias field, frot- Each initial condition contributes a particular amount to the variance in 
the phase between the relevant internal states. These contributions can be converted to partial contributions to a 
decoherence rate, and a numerical integral over a thermal distribution of collisional initial conditions can yield the 
total decoherence rate. We have pursued this program to a greater or lesser extent with the decoherence mechanisms 
discussed in this subsection and the one immediately following, but the results are not especially illuminating and we 
have used them primarily to confirm that the power-law expressions discussed below represent only overestimates of 
the decoherence rate, and that for experimental parameters of interest, the decoherence rate will be conservatively 
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less than 1 s"'^. 

The main question is whether T is high enough to include significant phase space for collision trajectories for which 
the peak value of fjon > ^rot (which is to say, large enough to transiently tip the direction of the total field by 
more than a radian). If so, then a single collision can cause decoherence and one can get a simple estimate of the 
cross-section for decoherence simply from the size of the impact parameter that leads to those events. There is a 
significant probability for collisions with £{on ^ Smt when 

T > 18 X f - " 

which leads to a decoherence rate 



T:>18K[ , (86) 

\5 V/cm/ 



T 



T 



1/2 



If T is instead so low that the Coulomb barrier suppresses collisions that could lead to a sufficiently large value of 
A and cause decoherence with a single collision, then decoherence will arise only from the combined effects of many 
collisions each causing small phase shifts that eventually random walk the science transition into decoherence. In this 
regime, the decoherence rate falls off very fast at low temperatures. For 



T< 18 X 



r \ 1/2 



5 V/cm 

typical collisions have f ion ^ i^^rot and the decoherence rate is 



« 0.13 X 



( ^ ^ ( JL.] (^12^ \ ^ 3- (89) 

VlOOO cm-3/ \ 15 K / \ 5 V/cm/ ' ^ ' 



1000 cm-3 / V 15 K y V 5 V/cm 



Both Eqs. [87] and |89] represent conservative estimates of the decoherence rate, and for an intermediate range of 
temperature, the decoherence rate will be less than whichever estimate gives the smaller value (Fig. [T^ . 

C. m-Level Changing Collisions 

A second source of decoherence can arise from ion-ion collisions that induce transitions between internal levels of a 
molecule. The dominant inelastic channel will be transitions between mp levels induced by a sufficiently sudden tilt in 
the quantization axis defined by the instantaneous local electric field. There are two conditions for such a transition 
to occur: (i) the direction of the total field must change by nearly a radian or more, so that there is significant 
amplitude for, e.g., an m_F = -1-3/2 level in the unperturbed electric field to suddenly have non- negligible projection 
on an mp = -1-1/2 level in the collision-perturbed field, and (ii) the time rate of change of the electric field direction 
must be comparable to or larger than the energy splitting between an m^? = 3/2 level and its nearest mp = 1/2 level 
in the field frot- 

Note that the first requirement is the same as the requirement for picking up an appreciable single collision Berry's 
phase. However, not all collisions that satisfy the first requirement will satisfy the second requirement. In particular, 
if the relative velocity in a collision is too low, then the time rate of change of the electric field direction will be too 
slow to satisfy the second requirement. Thus, given that the first requirement is satisfied, then the second requirement 
will not be satisfied when 

r<5Kx-^2^, (90) 
5 V/cm 

In this limit, the second requirement is more stringent than the first requirement, which means that the rate of 
m-level changing collisions will be smaller than the rate of single-collision Berry's phase-induced decoherence. In the 
opposite limit, we expect the second requirement will be met whenever the first requirement is met, and thus we would 
expect that the two channels of decoherence, m-level changing and single-collision Berry's phase, will be comparable 
in magnitude. 

Looking at particular collision trajectories in more detail, we see that there are trajectories that can cause an m-level 
change but for which there is no contribution to Berry's phase because the electric field traces out a trajectory with no 
solid angle (for instance, if the classical impact parameter b is parallel to frot)- We also note that our formulation of the 
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FIG. 12: Inverse coherence times, r~^, due to geometric phases accumulated during ion-ion collisions as a function of (a) 
collision energy in temperature units and (b) £rot- Dotted lines are approximations given in Eqs. [87] and 1891 Solid lines 
are more involved estimates based on integrals over collision parameters, but are still based on approximations so as to be 
conservative. The ion density was taken to be n = 1000 cm~'^. 



requirement of sweep rate for m-level changing collisions neglects the fact that £ion will not only change the direction 
of the total electric field (fion + ^rot) but also in general will change its magnitude. For most impact parameters, 
the magnitude of the total electric field will increase, thus suppressing nonadiabatic effects. However, a narrow range 
of impact parameters exists where the magnitude of the total electric field decreases, thus enhancing nonadiabatic 
effects. However, the above scaling laws account for the majority of collisions. 

In the end, we are less interested in the actual rates than we are in putting conservative limits on decoherence rates. 
For instance, in calculating the curves in Fig. I12[ we pessimistically took a worst-case geometry, fiot -L £iom which 
gives an upper limit on the size of the effect. Thus we estimate that: 

• For r < 5 K X g v/cm total coUisional decoherence, including both m-level-changing and Berry's-inducing 
effects, will be less than or equal to the value given by solid curves in Fig. [T^l while 

• For T > 5 K X ^ f//* , the total coUisional decoherence will be no greater than twice as large as the rate indicated 

5 V/cm ' ^ ^ 

by those curves. 

VI. CONCLUSIONS: PRECISION AND ACCURACY 

Recall from Sec. Uthe three components to the sensitivity figure-of- merit: 
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A. Coherence Time 

Conclusions of Sec. IIVI and |Vl Taking into account only collisional decoherence, and all the questions associated 
with being in rotating fields and in trapping fields, we would anticipate a coherence time longer than one second. 
Black-body thermal excitation of the J=l rotational level will also be well over one second. Vibrational black-body 
excitation for the v=0 state is estimated at 6 s for HfF"*" in a 300 K envir onm ent. Thus the dominant limitation to 
coherence will likely be the radiative lifetime of the '^Ai state, estimated jl38j at 390 ms for HfF+, and still longer 
for ThF+, for which the '^Ai state is predicted to be still lower in energy. The largest uncertainty in the lifetime 
calculation is the uncertainty in the "^Ai -> decay energy, calculated to be 1600 cm^^ in HfF+. 

B. fcff 

£cff in HfF+ is calculated by Meyer and coworkers to be 30 GV/cm and by Titov et al. to be 24 GV/cm [lOll |. 
For ThF+, Meyer calculates 90 GV/cm 28]. The uncertainties in these numbers are hard to assess, but they are 
very likely accurate to better than a factor of two and, if ongoing spectroscopic studies provide experimental values 
of hyperfine and fine structure that confirm the ab initio values predicted by the St. Petersburg group, our confidence 
in the precision of calculated SoS will be much higher. 

C. Count Rate and Summary of Expected Precision 

We are producing HfF+ ions by photoionization in a relatively narrow range of quantum states, and can estimate 
yield per quantum level within the desired trapping volume at perhaps 100 ions per shot, but we have just begun to 
characterize the efficiency of the process and very little optimization has been done. Our design efficiency for reading 
out spin states of trapped ions via laser- induced fluorescence (LIF) is 4%, but that has not been verified yet. With 
a large uncertainty, then, we may detect about one ion per shot with four shots per second. Overall, precision in 
one hour could be about 10 mHz. For ten hours of data, we anticipate (very roughly) a raw precision at 5 x 10^^^ 
e cm in HfF"*", and 1.5 x 10^^^ e cm in ThF+. We are investigating several more efficient alternatives to LIF for spin 
readout, including in particular resonantly enhanced photodissociation or second photoionization. Even if we detect 
as many as four ions in a shot, single shot precision will be no better than 300 mHz, which sets a relaxed requirement 
for suppressing experimental shot-to-shot noise. 

D. Accuracy, Systematic Error 

We have not completed a systematic study of the consequences of all possible violations of reflection symmetry in 
the trapping fields, but work in this direction is ongoing. 
For now, we make the following three observations: 

i) For the field asymmetries we have analyzed to date, realistic estimates for the magnitude in as-constructed field 
imperfections lead to systematic errors on the order of a few 10~^^ e cm or less. While this is not yet as accurate as 
our ultimate ambitions, it would represent roughly a factor of thirty improvement on the existing best experimental 
limit. 

ii) Asymmetries analyzed to date lead to systematic errors whose signs reverse when the direction of rotation Wrot 
reverses. If we combine measurements made with clockwise and counterclockwise field rotation, the errors vanish. 
Ideally, we'd like to design sufficient accuracy into the experiment so that the chop in field rotation is not needed to 
achieve desired accuracy, but as a practical matter we will of course run the experiment both ways, averaging the 
results to get ultimate accuracy, and differencing them to diagnose experimental flaws. 

iii) Auxiliary measurements are envisioned to characterize and shim out flaws in the as-constructed trap. For 
instance, we plan to be able to shim the equilibrium position of the ion cloud up and down along the trap axis, and at 
each location measure the energy difference E^-Erf. Unlike Eq-E;,, E^-Erf is highly electric-field sensitive. The result 
will be a precise measurement of any spatial gradient in ^rot • 

iv) All systematic errors we have analyzed to date have strong dependencies on quantities such as Wrot, Smt, ^rot, 
and on the ion-cloud temperature and density, and the trap confining frequencies. A true signal from a nonzero value 
of de will be largely independent of all those quantities. We anticipate making a number of auxiliary measurements 
with the experimental parameters tuned far away from their optimal values to deliberately exaggerate the size of 
systematic errors and allow us thus to characterize their dependencies in less integration time than that required for 
ultimate sensitivity. Even so, and as is often the case in precision measurement experiments, sensitivity and accuracy 
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are coupled. To the extent we can measure de to high precision at many combinations of experimental parameters, 
we will better be able to detect and reject false signals. 

We believe the experiment as we have described it should have the capability to improve the limit on the electron's 
electric dipole moment to 10~^^ e cm. As of this writing, the largest contribution to the uncertainty in our ultimate 
capability has to do with unknown efficiencies of state preparation and read out. More specialized publications from 
our group addressing progress in these areas are forthcoming. 

E. Summary 

Until now, molecular ions have not been viable candidates for symmetry violation searches largely due to the fact 
that applying electromagnetic fields to manipulate the internal states of the molecule would also violently perturb 
the translational motion of the ions. In this work, we have proposed a technique to overcome this obstacle - namely 
applying an electric field that rotates at radio frequencies. The specifics of performing high-resolution electron spin 
resonance spectroscopy under these conditions were analyzed. In particular, we have shown that a significant advance 
towards detecting the permanent electric dipole moment of the electron can be made by probing the valence electrons 
in a ground or metastable '^Ai level of an ensemble of trapped diatomic molecular ions. 

Note added in proof: Since the submission of this work, a new experimental limit on the electric dipole moment of 
the electron has been achieved using YbF molecules: \de\ < 10.5 x 10^^^ e cm jl39l |. 

Appendix: Typical Experimental Parameter Values 

de£eS — 2Tr X 0.36 mHz, transition energy between mp = +3/2 and tof = —3/2 states in ThF+ if — 1.7 x 10^^^ e 
cm. 

rfmf — +1.50 a.u. 27r X 2 MIIz/(V/cm), electric dipole moment of HfF+ in the molecular rest frame. 

i^^rot = 5 V/cm, rotating electric field. 

cjj-ot = 27r X 100 kHz, frequency of rotating electric field. 

Erot = 1800 K, typical kinetic energy in rotational micromotion. 

Trot = 0.6 mm, radius of circular micromotion. 

rfmf^rot = 27r X 10 MHz. 

(3/2)7j;-=3/2rfmt^^rot = 27r X 5 MHz, Stark shift of mp = ±3/2 states of '^Ai levels in rotating electric field. 
Wef = 27r X 10 kHz, A-doublet splitting between opposite parity '^Ai J=l states. 
5f=3/2 = 0.03, magnetic g-factor in '^Ai mp = ±3/2 states. 
Brot — 70 ^G, rotating magnetic field. 

^gpfiBBrot — 2t: X 8 Hz, Zeeman splitting between mp = +3/2 and mp — —3/2 states due to Srot- 

Sgp=3/2/ 9f=3/2 ~ 3 X 10^*, fractional difference of magnetic g-factor for upper and lower levels, for parameters 

shown in the Appendix. 

A « 27r X 2 Hz, splitting at the avoided crossing between mp — +3/2 and mp = -3/2 levels, for parameters shown 
in the Appendix. 

B± = 25 mG, anticipated scale of transverse magnetic field. 

r = 0.5 cm, characteristic rms radius of trapped ion cloud. 

T = 15 K, characteristic temperature of trapped cloud. 

uji = 2tt X 1 kHz, typical trap confining frequency. 

£2 = 5 mV/cm, typical axial electric field applied for confinement. 

frf = 75 mV/cm, typical Paul trap electric field strength, at typical cloud radius. 

< £r{ > — 5 mV/cm, typical radial confining electric field, averaged over one Paul cycle. 

cjrf = 27r X 15 kHz, typical "rf freq" for Paul trap. 

Ei-t = 15 K, typical kinetic energy in Paul micromotion. 

Ehf = 27r X 45 MHz, hyperfine splitting between F = 1/2 and F — 3/2 states of ^^Ai J=l level. 
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